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Abstract 

The discovery of magnetic and compositional effects in the low temperature properties of multi- 
component glasses has prompted the need to extend the standard two-level systems (2LSs) tunneling 
model. A possible extension [T] assumes that a subset of tunneling quasi-particles is moving in a three- 
welled potential (TWP) associated with the ubiquitous inhomogeneities of the disordered atomic structure 
of the glass. We show that within an alternative, cellular description of the intermediate-range atomic 
structure of glasses the tunneling TWP can be fully justified. We then review how the experimentally 
discovered magnetic effects can be explained within the approach where only localized atomistic tunneling 
2LSs and quasi-particles tunneling in TWPs are allowed. We discuss the origin of the magnetic effects 
in the heat capacity, dielectric constant (real and imaginary! parts), polarization echo and SQUID mag¬ 
netization in several glassy systems. We conclude by commenting on a strategy to reveal the mentioned 
tunneling states (2LSs and TWPs) by means of atomistic computer simulations and discuss the micro¬ 
scopic nature of the tunneling states in the context of the potential energy landscape of glass-forming 
systems. 
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1 Introduction 


The physics of glasses, especially at low temperatures, continues to attract considerable interest. At low 
temperatures glasses are believed to be characterized by low-energy excitations going under the name of 
tunneling systems, or states (TSs) which are normally described in terms of double-welled potentials (DWPs) 
and two-level systems (2LSs) with energy asymmetry and tunneling barrier uniformly distributed in the 
amorphous solid mm- Little is still known about the nature of the TSs, but the general consensus is still 
that the intermediate-range atomic structure of glasses is well described by Zachariasen’s 1932 continuous 
random network model [31[5] (thus, homogeneously disordered like for a liquid) and the 2LSs arise out of two 
slightly similar, localized atomic configurations. With this characterization the 2LSs have been employed in 
the 1970s and 80s to explain with some success the anomalies in the properties of glasses at low temperatures. 
Recently, the field has been witnessing a renaissance since the discovery that there are unexpected magnetic 
effects in insulating glasses [gizlIHlE] and that the TSs are responsible for the decoherence of phase qubits 
made of superconducting Josephson-junctions [M El m Haim. The 2LSs are thought to be ubiquitous in 
the junction’s tunneling barrier which, though very thin, is considered to be amorphous and described by 
Zachariasen-Warren’s model mm- On the other hand, the low-temperature glasses have also been studied 
as a paradigm where to conduct research on the physics of aging so common to many disordered systems out 
of equilibrium [ElEl]- this context, the relevant degrees of freedom have also been described in terms of 
2LSs, with some success [T7] . 

Yet, deviations from the behaviour predicted by the standard tunneling model (STM) have challenged the 
validity of the model in the case of multi-component glasses with variable content of the good crystal¬ 
forming (GCF) component (e.g. (Si02)i-x(K20)i_a; with changing x [ElilS]) and especially in glasses of 
compositions like Ba 0 -Al 203 -Si 02 in the presence of a weak magnetic field [20]. In such glasses (the mixed 
alkali-silicates have not yet been investigated in a field, but we predict important magnetic effects there, 
and x-dependent, too) a puzzling non-monotonous magnetic-field dependence has been revealed in most 
physical properties 0 El [i 0. The magnetic effect is normally weak, but orders of magnitude larger than 
expected from basic thermodynamic considerations. The STM is unable to explain the compositional and 
magnetic effects, thus a suitable extension of the tunneling model for both situations has been proposed by 
one of us [1]. This extended tunneling model (ETM) rests upon the existence in the multi-component (or 
contaminated mono-component) glasses of regions of enhanced regularity (RERs) in the atomic structure 
of the incompletely-frozen (in fact) amorphous solid. A complete mathematical description and physical 
justification of the ETM thus calls at the very least for a partial demise of the Zachariasen-Warren’s vision 
of the intermediate structure of glasses. 

Research on the mechanism of the glass phase transition is also witnessing intense activity |21j , with almost all 
of the investigations conducted from the high-temperature (liquid or supercooled-liquid) side of the transition. 
Any attempt at a new look at (and demise of) the Zachariasen-Warren’s continuous random network model 
description of the atomic structure of glasses from the low-temperature side is therefore likely to provide new 
insight also on the nature of the glass transition itself. 

In this paper we provide a new scenario for the atomic structure of glasses, the cellular model, a description 
within which the phenomenological assumptions for the ETM’s mathematical formulation become completely 
justified. The cellular model provides for a more realistic mathematical formulation in terms of a tetrahedric 
four-welled tunneling potential, the triangular three-welled version employed being nothing but a poor-man’s, 
but probably very realistic, version affording a much simplified mathematical description. Within the cellular 
approch to the structure of glasses the most significant local tunneling potentials (also in terms of number 
density) turn out to be the DWPs, for a single (or very few) atomic particles, and the tetrahedric four-welled 
potential (TFWP) for a correlated cluster of S> 1 charged atomic particles. A reasonable and very useful 
simplification for the TFWP is then the replacement of the N interacting and tunneling atomic particles with 
a single quasi-particle subject to a triangular tunneling potential (TWP) carrying renormalized parameters 
(charge, magnetic threaded area, energy asymmetry and tunneling probability) moving about one face of the 
tetrahedric potential. The renormalization is fully justified by the proximity of four similarly quasi-ordered, 
close-packed atomic cells (RERs) and the reasonable assumption that most of the charged particles will avoid 
the interstice’s centre. We present therefore the mathematical description of the TWP and its quantum 
mechanics in the appropriate limits for practical applications. We apply the results for the description of 
the density of states and the temperature and magnetic-field dependence for the specific heat of the glass, 
then for the dielectric constant (real and imaginary part) in the linear regime, and for the polarization echo - 
always in the presence of a magnetic field. Finally we apply the model to the description of the temperature 
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and magnetic-field dependence for the SQUID-magnetization, showing that the ETM keeps providing a good 
description for this quantity up to room temperatures. As a check that the tunneling states we use are 
realistic representations of the remnant degrees of freedom in the disordered solid phase we also outline a 
strategy for the study of these local tunneling states in the potential energy landscape (PEL) with the help 
of classical numerical simulation of a large number of particles interacting with a BLJ potential. 

The paper is organized as follows. In Section 2 we introduce the cellular model for the structure of glasses, an 
important refinement of Zachariasen’s continuous random-network model, and examine the likely tunneling 
states that emerge from such cellular picture, concluding that only DWPs and TFWPs should be relevant 
for the physics of glasses below the glass transition temperature Tg. In Section 3 we review the relevant 
quantum mechanics of the three-welled, poor man’s version of the TFWP, a version that has been used so 
far to obtain a reasonable single explanation for all the anomalies (and deviations from STM predictions) 
due to composition changes and to the magnetic field. We also show how to evaluate the magnetic density 
of states (DOS) n{E,B) and, in Section 4, we examine the magnetic-field dependent heat capacity Cp(T,B) 
comparing with some published data for the multi-silicates. In Section 5 we do the same for the dielectric 
constant, real part e' and imaginary part e", also showing some comparison with available data at low (kHz) 
frequency. In Section 6 we examine the application of the ETM to the explanation of some of the data for the 
polarization echoes in the silicates and in glycerol, also explaining the so-called isotope effect. In Section 7 we 
examine the question of the paramagnetism of insulating glasses per se, that is that due to the tunneling states 
themselves and which we find to be stronger than the Langevin paramagnetism of the trace paramagnetic 
impurities contained in the multi-silicates. Section 8 contains a discussion about guiding principles on how 
one should proceed, within a molecular dynamics simulation of a model glass former (the BLJ system in 3D, 
for example), in order to characterize the relevant tunneling states from a numerical analysis of the PEL and 
thus make contact with the original phenomenological modelization at low temperatures. Finally, Section 9 
contains a Summary of the results obtained and our Conclusions. 


2 The Cellular Model for the Atomic Structure of Glasses and the 
Three-Welled Tunneling Potential 

What is a glass? According to Shelby “A glass can be defined as an amorphous solid completely lacking 
in long range, periodic atomic structure and exhibiting a region (of temperature) of glass transformation 
behaviour” [^. This, however, does not mean that the atomic arrangement of a glass should be the same 
as that for a liquid, as is implicit in Zachariasen’s 1932 proposed continuous network picture mis], which 
has been widely adopted by scholars (at least in the West, see below). This picture differs from that of a 
liquid only in that a dynamical arrest has occurred, without specifying its ultimate origin. In a spin-glass 
the ultimate origin of dynamical arrest is magnetic frustration (with or without disorder), but for ordinary 
structural glasses it remains mysterious and an important open issue EH na mils]. Prior to Zachariasen’s 
scheme, however, the Soviet scientist A.A. Lebedev had proposed, in 1921, the concept of “crystallites” [26] . 
small undefined crystal-like regions clumpt against each other in random orientation to contain altogether 
all atoms in the substance. Later, Randall proposed that these be real micro-crystals and explained the 
rounded-up X-ray spectra from glasses in this way m- However, the density of glasses is some 10% less than 
that of poly-crystalline aggregates and the thermal properties of glasses also cannot be explained through the 
Lebedev-Randall picture. Despite these observations, the West-Soviet controversy continues to these days. 
The Zachariasen-Warren model of glass structure was in fact criticized by Hagg [28] in the West right in the 
early days of X-ray crystallography and a good review of the status quo of this controversy has been recently 
provided by Wright [29] who concludes from a re-analysis of X-ray and neutron-scattering data from many 
covalent-bonded and network glasses that indeed cybotactic groupings (better-ordered regions) may well be 
present and frozen-in in most glasses, especially if multi-component. The formation of polyclusters, instead 
of crystallites, in most glasses is the latest claim by the Soviet school [SO] |3T], based on observations and 
thermodynamic reasoning. As nicely set out by Bakai [30] . the incipient crystals forming at and below 
(melting point) are in constant competition with kinetically swifter (for glass-forming liquids) polyclusters 
(or crystallites) that can win thermodynamically and kinetically over crystals during a rapid enough quench. 
For e.g. the oxyde glasses (e.g. window glass, of composition Ca0-Na20-Si02) the dilemma is most simply 
put in these terms. Pure good glass-former (GGF) Si02 melts at around 1900 K and vitrifies at around 1475 
K, whilst pure GCFs CaO and Na20 melt at 2886 K and 1405 K, respectively. Yet, the mixture GaO-Na20- 
Si02 vitrifies at around 750-850 K (depending on % composition) and that is way below the melting point of 
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very bad glass-formers CaO and Na20. Since phase separation is commonplace in such ternary mixtures (the 
recent X-ray tomography imaging of the related molten system Ba0-Na20-Si02 is compelling [31]), what is 
then stopping the pockets and channels formed by CaO and Na20 within the GGF Si02’s network to at least 
initiate some degree of crystalline ordering? In the mono-composed oxydes it will be up to the impurities to 
take the place of the GGF phase separation. 

On the experimental side, the concept of de-vitrification is gradually taking shape with reports of metallic 
glasses [33] and also monocomponent high coordination covalent solids like Si forming paracrystals in the 
amorphous solid |34j . Therefore, the stance will be taken in this paper that only the purest mono-component 
glasses may abide to the Zachariasen-Warren continuously random-network model of a glass, whilst the vast 
majority of real glasses will be organized otherwise at the intermediate-range atomic structure. 

The structure proposed is a cellular-type arrangement of better-ordered regions (regions of enhanced order¬ 
ing, RERs) that can have complicated, maybe fractal, but compact shapes with a narrow size distribution 
and interstitial regions between them populated by still fast-moving particles (normally charged, probably 
dangling bonded ions). We remark that a cellular-type structure was already embodied by the “crystallite” 
idea of Lebedev, but now no micro-crystals are here claimed to exists (except perhaps in the ceramic-glasses 
case, to an extent). The RERs are more like Wright’s “cybotactic regions” [29| or Treacy’s “paracrystals” 
[34] and issue from the dynamical heterogeneities (DHs) picture of the glass-forming supercooled liquid phase 
above Tg (see e.g. Fig. [5 [lll|33|3a|33. The DHs picture recognises the presence of regions of “slower” 
and “faster” particles, and inspection of the slower-particle regions in the supercooled liquid reveals that 
these are also better ordered (solid-like) whilst the faster-particle regions are much more liquid-like. DHs are 
ubiquitous in most supercooled liquids [37] and the claim (yet still somewhat speculative) here is that the 
slower-particle regions will grow on approaching Tg, but only up to a finite size and will be giving rise to the 
RERs in the frozen, glassy phase. Simulations in the frozen glassy state of the slower- and faster-particle 
regions confirm that a DHs picture applies also below Tg |38j and with the slower-regions increasing in size 
as T —)■ 0. Earlier simulations (always for model systems) [33] pointed out the difficulty of simulating the 
DHs picture below Tg and came up with a picture of the slower-regions growing, possibly diverging in size, 
on approaching Tg. In this paper, however, the stance will be taken that the slower-regions’ average size 
growths, in real systems, but does not diverge at Tg or at any other characteristic temperature. A full 
numerical proof of this fact is, however, still lacking. We remark that a cellular structure for the glasses 



(a) (b) 

Eigure 1: (colour online) (a) Map of the relaxation times of one MD simulation run of a 2D system of 780 soft 
disks, interacting via a repulsive 1/r^^ potential cut off at some distance (from Ref. [3S]). Circles represent 
particle positions, the slowest 40% coloured black, the fastest 40% shaded gray and the intermediate 20% 
are unfilled circles. Note the distinct clustering of slow particles (dynamical heterogeneities, DHs). (b) The 
slower regions have been schematically highlighted to show their incipient cellular structure. 

had been proposed in the past by de Gennes |40j and, in the context of the low-temperature anomalies, by 
Baltes m who was able to explain the linear in T anomaly in the heat capacity Cp (but not those in the 
acoustic properties which require the introduction of 2LSs). A very similar picture is that of the polyclusters 
of Bakai [3D], as already mentioned, but here the thermal history of the cells, or RERs, is ascribed directly to 
the DHs situation already present above Tg. In this approach the RERs, like grains in the frozen structure 
below Tg, are the thermal history continuation of the slower particle regions of the DHs above Tg and contain 
the bulk of the 2LSs. These 2LSs are atomic tunneling states arising from the grains’ own disorder causing 
a distortion in the natural bond angle between neighbouring cations with the anions becoming susceptible 
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of taking up two nearly equivalent positions differing by a tiny energy (see Section 8). Some 2LSs may be 
located at the meeting point between two grains, or cells (two RERs). In the interstitials between grains the 
faster particles of the DHs eX T > Tg give rise for T < Tg to regions where a large number N (on average) of 
charged atomistic tunneling particles are costrained to move in a coherent fashion due to the high Coulomb 
repulsion between them. Fig. shows in a schematic way how the atomic/ionic matter can get organised 
below Tg in a real glass. Since the charged ions (dangling bonds, most likely) should act as a coherently 



Figure 2: (colour online) (a) 2D cartoon of the cellular structure of an amorphous solid just below Tg. The 
RERs (black-circled blobs, an oversimplified schematics for fractal-like, but compact objects) have grown to 
completely fill the space and enclose atomistic tunneling states of the 2LS type (blue blobs). At the same 
time, in the RER interstitials (yellow regions, connecting to each other) the trapped, charged and faster 
particles of the DHs existing above Tg and now probably charged dangling bonds, give rise to coherently 
tunneling large groups of ions to be represented by a single fictitious quasi-particle (orange dot) subjected to 
an effective tunneling potential having four natural wells in distorted tetrahedral configuration (b) for close- 
packed RERs. (b) The tetrahedral four-welled potential (TFWP) in a 3D representation with colour-coded 
potential intensity (dark=deepest, light=highest). 


tunneling ensemble, it seems natural to simplify the description of the physics at the lower temperatures 
using phonons, propagating in the collection of cells now jammed against each other, and remnant localized 
degrees of freedom acting as TSs. These TSs will be the 2LSs within the cells and at their points of contact 
(owing to inherent disorder in the cells’ atomic arrangement) and effective quasi-particles sitting in the close- 
packed cells’ interstices and representing the collective motion of the coherently- tunneling ions trapped in 
each interstice (Fig. [^. The quasi-particle will be subjected to an effective potential of distorted- tetrahedral 



Figure 3: (colour online) (a) Contour plot of a possible realization of the 2D effective three-welled potential 
(TWP) most likely felt by the quasi-particle of those charged real particles dangling from a group of three 
RERs on each one face of the tetrahedral configuration of an interstice formed by close-packed RERs well 
below Tg. (b) 3D visualization of the same example of a TWP potential. 


shape characterized by four wells for each interstice, with a high barrier in the interstice’s centre. De facto 
this 3D interstitial TFWP potential can be replaced with four local 2D potentials for the four quasi-particles 
describing the coherent tunneling of the dangling-bond particles sitting near each face of the tetrahedron. 
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close to a group of three (on average) quasi-ordered cells (Fig. [^. Because of the better ordering implicit in 
this model of the glassy intermediate-range atomic structure and in each cell, the three wells of each effective 
local 2D potential for the tunneling quasi-particles (four per interstice, on average) will be near-degenerate in 
terms of their ground-state energy asymmetries: Ei ~ E 2 — ^ 0. With this qualitative picture in mind, 

we now turn to the mathematical description of the physics of the remnant, still ergodic degrees of freedom 
(phonons in the cells’ network, 2LSs and ATSs (anomalous tunneling systems, four in each interstice). For 
our model of a real glass, by construction the 2LSs will be more numerous than the ATSs. 


3 The Quantum Mechanics of the Three Welled Potential 

In this approach [1] the relevant degrees of freedom, beside the phonons, are dilute collections of independent 
2LSs, described by the STM, and of fictitious quasi-particles tunneling in TWPs. 

The formulation of the STM (the 2LS model) for the low temperature properties of glasses is well known. 
One assumes a collection of DWPs distributed in the substance and represented each by a 2x2 Hamiltonian 
of the form, in the potential-well (or real-space) representation: 


n^Ls = -2 


A Ao \ 
Ao -A J ■ 


( 1 ) 


Here the parameters A (the energy asymmetry) and Ag (twice the tunneling parameter) are typically charac¬ 
terized by a probability distribution that views A and In(Ao) (the latter linked to the DWP energy barrier) 
broadly (in fact uniformly) distributed throughout the disordered solid [i2] : 

'P 2 Ls(A, Ag) = — (2) 

where some cutoffs are introduced when needed and where P is a material-dependent parameter, like the 
cutoffs. In fact Eq. § embodies the Zachariasen-Warren hypothesis for the intermediate atomic structure 
of a glass, assuming broadly distributed the energy asymmetry A = E^'^ — E^^\ hence the single-well ground 
state energies E^'^ and E^^'^ themselves, as well as the potential barrier height Vg appearing in relations such 
as: 
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where the first relation is the generic WKB result for an arbitrarily shaped DWP (where m is the particle’s 
mass, D its single-well harmonic frequency (tunneling attempt frequency) and d the tunneling distance) and 
the second formula refers to a symmetric (A = 0) DWP made up by two superimposed parabolic wells. In 
fact, the distribution ([^ in the end refers to the combination of parameters Ag/SD. The energies of the two 
levels |0 > and |1 > are then £g_i = ±4 a/A^~+~A§ and so on [21142]. 

The tunneling Hamiltonian of a particle in a TWP is easily written down, in the same low-T spirit as for a 
2LS, as a generalization of the above matrix formulation to three levels [Tj: 


T~1-3LS = 
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where Ei,E 2 ,E^ are the energy asymmetries between the wells and Dq is the most relevant tunneling 
amplitude (through saddles of the PEL, in fact). This 3LS Hamiltonian has the advantage of readily allowing 
for the inclusion of a magnetic field B > 0, when coupling orbitally with a tunneling “particle” having charge 
q {q being some multiple of the electron’s charge —e) [T|: 


TdsLsiB) = 


E 2 DoE'^/^ 

Doe-^v/3 Es ) 
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( 6 ) 








where tf/i is the Peierls phase for the tunneling particle through a saddle in the field, and (f is the Aharonov- 
Bohm (A-B) phase for a tunneling loop and is given by the usual formula: 


„ $ 

f = 2 x 5 ;, 




(7) 


$0 being the appropriate flux quantum {h is Planck’s constant) and d) = B-Sa the magnetic flux threading the 
area S'a formed by the tunneling paths of the particle in this simple (poor man’s, yet as we have seen realistic) 
model. The energy asymmetries Ei, E 2 , E 3 typically enter through their combination D = ^yE1 + E 2 + E'^. 
For 71^=3 wells an exact solution for the k=0, 1, 2 eigenvalues of the multi-welled tunneling Hamiltonian 
Eq. (§ is still possible: 


£k = 2 D 0 1 /1 ~ 
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9k = 0,-|-|7r, — Itt distinguishing the three eigenstates. In the physically relevant limit, which we now 
discuss, in which 0 (weak fields) and D = ^E^ + E 2 + E'^ —>• 0 (near-degenerate distribution), and at 
low temperatures, we can approximate (in a now simplified calculation) the riu, = 3 - eigenstate system with 
an effective 2LS having a gap AS = Si — Sq widening with increasing ip it Dq > 0 (see below): 


limA£ ~ 

73 


\j+ 2^EI + E 2 + El) 


£>2 


(9) 


(a trivial rescaling of Dq and of the Ei was applied). One can easily convince oneself that if such a TWP is 
used with the standard parameter distribution, Eq. ([^ with D, Dq replacing A, Aq, for the description of the 
TS, one would then obtain essentially the same physics as for the STM’s 2LS-description. In other words, 
there is no need to complicate the minimal 2LS-description in order to study glasses at low temperatures, 
unless structural inhomogeneities of the RER-type and a magnetic field are present. Without the RERs, 


hence no distribution of the type (10) below, the interference from separate tunneling paths is only likely to 


give rise to an exceedingly weak A-B effect. Hence, it will be those TSs nesting between the RERs that will 
give rise to an enhanced A-B interference and these TSs can be minimally described - most appropriately - 
through Hamiltonian Q and with distribution of asymmetries thus modified to favour near-degeneracy (P* 
being a dimensionless material parameter) [1]: 


V*qls{Ei,E2,Eq-Dq) = 


p* 


Dq{EI + El + El) ■ 


( 10 ) 


We remark that the incipient “crystallinity” of the RERs calls for near-degeneracy in Ei,E 2 ,Eq simultane¬ 


ously and not in a single one of them, whence the correlated form of (10). We now have three-level systems 
(3LSs) with energy levels Sq < Si ^ £ 2 , periodic in p. The typical spectrum, with Dq > 0 (see below), is 
shown in Fig. |^as a function of p and one can see that the third, highest level S 2 can be safely neglected for 
most applications. Other descriptions, with TFWPs or modified three-dimensional DWPs are possible for 
the TSs nested in the RERs and lead to the same physics as from Eqs. (1^ and (101 above [33] (which describe 


what we call the anomalous tunneling systems, or ATSs, nesting in the interstitials between the RERs). The 
final and most important consideration for the construction of a suitable mathematical model is that the 
TSs appear to be rather diluted defects in the glass (indeed their concentration is of the order of magnitude 
of that for trace paramagnetic impurities, as we shall see), hence the tunneling “particles” are embedded 
in a medium otherwise characterized only by simple acoustic-phonon degrees of freedom. This embedding, 
however, means that the rest of the material takes a part in the making of the tunneling potential for the TS’s 
“particle”, which itself is not moving quantum-mechanically in a vacuum. Sussmann |44j has shown that this 
leads to local trapping potentials that (for the case of triangular and tetrahedral perfect symmetry) must be 
characterized by a degenerate ground state. This means that, as a consequence of this TS embedding, our 
poor man’s model, Eq. (§, for the ATSs must be chosen with a positive tunneling parameter [I]: 


Dq > 0 


( 11 ) 
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Figure 4: The energy spectrum {Dq = 1 units) of the ATS (3LS, TWP) model, in the appropriate limits of 
weak field and near-degeneracy due to the embedding within RER’s interstitials. On the horizontal axis the 
A-B phase (p (x B. 
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Figure 5: (colour online) The magnetic-sensitive part of the density of states (DOS) as a function of the 
energy gap E and different A-B phases p (proportional to the magnetic held B) {uatsP* has been set to 
1). The rapid shift of quantum states to higher energy when a very weak B is switched on is the physical 
explanation for the origin of the magnetic effects. 


where of course perfect degeneracy is always removed by weak disorder in the asymmetries. The intrinsic 
near-degeneracy of (10) implies that this model should be used in its D/Dq ^ 1 limit, which in turn 
reduces the ATSs to effective magnetic-held dependent 2LSs and greatly simplihes the analysis together with 
the limit > 0 which we always take for relatively weak magnetic fields. The ETM, first proposed in 
[T], consists then in a collection of independent, non-interacting 2LSs described by the STM and also 3LS 
TWPs, described by Eqs. ^ and (10) above, in the said D/Dq 1 and p ^ 0 limits, the 3LSs nested 
in the interstitials between the close-packed RERs and the magnetic-held insensitive 2LSs distributed in 
the remaining homogeneously-disordered matrix of RERs and at their touch points or interfaces |46j . In 
Fig. we illustrate the behaviour of the density of states (DOS) for this model as a function of the gap 
energy E for different p values. This hgure shows the physical origin of the magnetic effects: the quantum 
states being conserved in number, they are very rapidly shifted towards high values of the energy when a 
magnetic held, even very weak, is switched on. Our ETM has been able to explain the magnetic effects in the 
heat capacity in the real [45] and imaginary m parts of the dielectric constant and in the polarization 
echo amplitude [47] measurements reported to date for various glasses at low temperatures, as well as the 
composition-dependent anomalies nms]. The new physics is provided by the magnetic-held dependent TS 
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Figure 6: (colour online) The magnetic-sensitive part of the density of states (DOS) as a function of the 
A-B phase ip (proportional to the magnetic field B) and different energies {uatsP* has been set to 1). 
The shape of this part of the DOS (coming from the TWPs with a parameter distribution (10) favouring 
near-degeneracy) is the ultimate source of all the magnetic effects. The cusp is an artifact of the effective 
2LS approximation [1], but also of the existence of upper and lower bounds for Dq owing to the nature of 
the RER glassy atomic structure. 


DOS, which acquires a term due to the near-degenerate TWPs [T] that gets added up to the constant DOS 
from the STM 2LSs (having density n 2 Ls)- 


gtot{E,B) = n 2 LsP + nATS-j^fATs{E, B)9{E — Ed) 

= 92Ls{E) + gATs{E,B) (12) 


where uats is the ATSs’ concentration, /ats is a magnetic-field dependent dimensionless function, already 
described in previous papers [T], and Ed is a material and B-dependent cutoff: 


9ats{E,p)= J U^dE^SiEjEj) J dDoP^LsiiPk}, Do)SiE - A£] 

' ^ 
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E < Ed- 


(13) 


Here, after a suitable renormalization of parameters: Ed = E^ = + E'omax'fi^y 

Dmin, Domin and D^max being suitable cutoffs (material parameters). The 1/E dependence of the ATS DOS 
is a consequence of the chosen tunneling parameter distribution, Eq. (10), and gives rise to a peak in gtot 
near Ec 2 that is rapidly eroded away as soon as a weak magnetic field is switched on. The form and evolution 
of the magnetic part of the DOS is shown in Fig. for some typical parameters, as a function of oc i? 
for different values of E. This behaviour of the DOS with B is, essentially, the underlying mechanism for 
all of the experimentally observed magnetic held effects in the cold glasses within this model: the measured 
physical properties are convolutions of this DOS (with appropriate H-independent functions) and in turn 
reproduce its shape as functions of B. 


4 The Magnetic Field Dependent Heat Capacity 

4.1 Theory 

As a hrst example, the total TS heat capacity is given by 

CpTs{T,B)= dE gtot{E,B)CpoiE,T) (14) 

Jo 
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where 


WT) = *,(^)%„h-(^) (15) 

is the heat capacity contribution from a single TS having energy gap E and where gtot{E, B) is given by Eq. 

In this Section we re-analyze some of the available data [35] for the magnetic effect in the heat capacity of 
two multi-component glasses, commercial borosilicate Duran and barium-allumo-silicate (AlBaSiO, or BAS 
in short) glass, in order to show that the ATS model works well for the magnetic-field dependent Cp. A 
systematic experimental study of Cp{T, B) around and below 1 K in these multi-silicate glasses was carried 
out by Siebert |35| and those data have been used, upon permission, by one of us [T| as the very first test of 
the above (Sections 2, 3) ETM [Tj. That earlier analysis best-fitted the Cp data by Siebert with the sum of 
Einstein’s jphT^ phonon term plus the 2LS 72 LsT non-magnetic contributions, as well as with Langevin’s 
paramagnetic and the ATS contributions (see below). The analysis came up with concentrations hj ~48 
ppm and, respectively, hj ~ 20 ppm instead of the quoted [48] 126 ppm (or 180 ppm in a different study |49|) 
and 102 ppm for Duran and for BAS glass, respectively. 

In order to better understand this large discrepancy we re-analysed Siebert’s data for Cp{T,B), after sub¬ 
traction of the data taken at the same temperatures for the same glass, but in the presence of the strongest 
applied magnetic field (8 T) |50j . In this way only the magnetic-field dependent contributions should remain 
in the data for Cp{T, B) = Cp{T, B) — Cp{T, oo). The parameters involved when fitting data are the cutoff 
Dmin and combinations of cutoffs, charge and area DominqS and DomaxqS [T], as well as: nFe^+ (Fe^"'' impu¬ 
rity concentration), n^g3+ (Fe^"'" impurity concentration) and uats (ATS concentration (always multiplied 
by P*)). The data from [35] (we have restricted the best fit to the three temperatures having the most 
data points around the peak of Cp{B)) have been best-fitted by using the following magnetic-dependent 
contributions: 


a. the known Langevin contribution of the paramagnetic Fe impurities (Fe^“'" and Fe^“'") having concen¬ 
tration nj'. 


Cj[T,B) = nj 


ksz^ 


1 


sinh ; 


2J+1 


sinh 


{2J+l)z 


(16) 


where z = and where g is Lande’s factor for the paramagnetic ion in that medium, is Bohr’s 

magneton and J the ion’s total angular momentum (in units fi.=l); kg is Boltzmann’s constant. We 
have assumed the same values the parameters g and J take for Fe^“*' and Fe^’*' in crystalline Si02: J=2 
with g=2 and, respectively, J = 5/2 with g=2 (we have adopted, in other words, complete quenching 
of the orbital angular momentum m, consistent with other Authors’ analyses 


b. the averaged contribution of the ATSs [T], written in terms of a sum of individual contributions from 
each ATS of lowest energy gap E (making use of Eqs. (121, (13) and 


Cats{T,v^) = -^-^ 




dE 


E 


cosher ^ 
cosn \2kBT 

E 


In 




/ dE - F 
'Ea 2 cosh ( 2 kBT 


In 


Domax \ E^ — 


(17) 


) ^\DominJ E'^ D^max^^ 


1 } 


or, re-written in a dimensionless form as 

Cats{T, v) = Cq{T, tp) 2T:P*nATskg^ [l{x cl) -^(^c2)] 

; ^mm) ^(^c2j^mm) “I” ^(^cl; ^Omm p) T(Xc 2? :^0max7’) J" 


( 18 ) 


where: 

• Ed 


2 

min 


= Vd. 

_ £^c1,2 

Xcl,2 — 2/CijT’ 


+ and Fg2 = 




Dry 


2 /cbT ’ 


etc.; 
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BAS glass 

Concentration [g 

Concentration [ppm] 

Tlp^2-\- 

i.oexiQi^ 

14.23 

Tlp^3+ 

5.00x10^® 

6.69 

P*nATS 

5.19x101® 

- 


Table 1: Extracted parameters (from the heat capacity data) for the concentrations of ATSs and Fe-impurities 
for the BAS glass. 


Temperature [K] 

^min [K] 

Dormn\i\S [Kl^j 

Doraa.\7\S [Kl^j 

0.60 

0.49 

4.77x10^ 

3.09x10® 

0.90 

0.53 

5.07x10^ 

2.90x10® 

1.36 

0.55 

5.95x10^ 

2.61x10® 


Table 2: Extracted tunneling parameters (from the Cp data) for the BAS glass. 


• /(a;) = xtanhx — Incoshcc; 

• ^(^>«) = Xrd2/^^ln(2/2-a2). 

and the following expression: 

Co(r, (p) = 2nP*nATskB In ({ 1^(2) - I{x,2).} (19) 

The angular average over the ATS orientations is performed by replacing ^ (in other words 

averaging cos^ 9, 9 being the orientation of S with respect to B). 

The above formula for Cats is actually correct only for weak magnetic fields (up to about IT). For higher 
magnetic fields one must make use of an improved form |53j for the ATS lower energy gap. In practice, this 
consists in making the following replacement for the A-B phase ip 

where the second expression holds after orientational averaging and where B* is the upturn value of the 
magnetic field for, e.g., the B-dependence of the dielectric constant e'(T, B) [55] . 


4.2 Comparison with available data 

The concentrations of the ATSs and Fe-impurities extracted from the best fit of the heat capacity as a function 
of B, for the BAS glass, are reported in Table[^ having fixed the concentrations, it was possible to extract the 
other parameters for the BAS glass (Table [^. The best fit of the chosen data is reported in Fig. [^a). The 
concentrations of the ATSs and Fe-impurities extracted from the best fit of the heat capacity as a function 
of B, for Duran, are reported in Table[^ having fixed the concentrations, it was possible to extract the other 
parameters for Duran (Table |^. The fit of the chosen data is reported in Fig. Bb). The first comment 
that we make is that these good fits, with a smaller set of fitting parameters, repropose concentrations and 
tunneling parameters very much in agreement with those previously obtained by one of us [T]. The problem 
with the concentrations of the Fe-impurities reported in the literature is that they do not allow for a good 
fit of the Cp{B) = Cp{B) — Cp{oo) data in the small field region when the Langevin contribution alone is 
employed (Eq. without Eq. (HI))- See for instance Section 5, Fig. The Langevin contribution drops 
to zero below the peak, whilst Siebert’s data definitely point to a non-zero value of C'p(O) = C'p(O) — C'p(oo) 
at B = 0 for any T > 0. This non-zero difference is well accounted for by our ETM and it comes from the 
ATS contribution to the DOS in Eq. (12). The parameter B* was left indetermined in these fits, owing to 


the fact that the measured Cp becomes very small at higher fields. 

The results of our Cp analysis definitely indicate that the concentration of paramagnetic impurities in the 
multi-silicate glasses is much lower than previously thought and extracted from SQUID-magnetometry mea¬ 
surements of the magnetization M{T, B) at moderate to strong field values and as a function of T. Therefore 
we will turn our attention to a re-analysis of the SQUID-magnetometry data in Section 7. 
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xio’® Duran 



(a) 


(b) 


Figure 7: (colour online) The heat capacity best fit for the a) BAS (or AlBaSiO) and b) Duran glasses. The 
continuous lines are the predictions from the ETM. 


Duran 

Concentration [g 

Concentration [ppm] 

Tlp^2-\- 

Tlp^3+ 

P*nATS 

3.21x101^ 

2.11x101^ 

8.88x10^® 

33.01 

21.63 


Table 3: Extracted parameters (from the heat capacity data) for the concentration of ATSs and Ee-impurities 
for Duran. 


Temperature [K] 

Dmin [K] 

Do^^n\i\S [KA^] 

Doma.\i\S [Kl^j 

1.11 

0.34 

4.99x10^ 

2.68x10*^ 

1.23 

0.32 

5.30x10-* 

2.50x10^* 

1.36 

0.32 

5.54x10-* 

2.46x10'* 


Table 4: Extracted tunneling parameters (from the Cp data) for Duran. 
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5 The Magnetic Field Dependent Dielectric Constant 


5.1 Dielectric properties of cold glasses: general theory in zero magnetic field 

To measure the frequency-dependent dielectric properties of the cold glasses one applies an ac electric field 
to the sample, typically at radio frequencies (RF). The linear response regime only will be considered in this 
work. The TSs then couple to this field via the electric charge or dipole moment of the tunneling “particle”. 
The applied electric field both modulates the energy splitting of the tunneling states and excites them from 
thermodynamic equilibrium. The electric field only affects the asymmetry energy A [2] . The influence of the 
electric field on the tunnel splitting Aq is usually neglected . The coupling of a 2LS to the external field 
therefore causes resonant processes like resonant absorption and stimulated emission. In the presence of the 
external electric field F the Hamiltonian matrix takes the form (coordinate representation): 


i? = i/o + Po • Fcr^ 


1 / A — 2po • F Aq 

2 V Ao -A + 2po • F 


( 21 ) 


Here po denotes the electric dipole moment of the fictitious particle, F = F^coswt is the time-dependent 
electric field. Diagonalizing the Hamiltonian (Eq. (H])) one can get (energy representation): 

= _i^) + (i t)po-Fcosa;t (22) 

The dynamics of the 2LSs is given by the change in the expectation values through the Bloch equations, 
which were first derived by Bloch in the context of magnetic resonance m- 

l{SyB,-S.By) 

'^=-Y^Sy+^[SA-SM (23) 

^ - ^ - (5.)) + 7 {S.By - SyB,) 

Here we have introduced the pseudo-spin f operator S= crj^^ where cr are Pauli’s matrices (in the 2LS 
energy representation), Ti is a characteristic time for the equilibration of the level populations of the 2LS, 
and T 2 is the transverse dephasing time due to spin-spin (i.e. 2LS-2LS) interactions. Also, {Sz) is the 
thermal equilibrium value of Sz given by {Sz) = \ tanh ^ ~ 7 is the appropriate (fictitious) 

gyromagnetic ratio and B= BdcS- Bac is a fictitious effective field made up of a static (dc) part and of an 
oscillating (ac) one, both proportional to the electric field with frequency lo. The dimensionless fictitious 
spin S processes around this fictitious effective field B, given by {fj, = hr): 


/iB 


2Ao 

~E~ 


Po • F, 0, E + 


2A 

1 ^ 


Po • F 


(24) 


Since the ac field is a small perturbation, one can expand (in linear response theory) (Sz) in a Taylor 
series by keeping terms up to the first order in Bac- The solution to the Bloch equations takes the form 
S{t) = S°(<) -I- S^(t), where S°(<) is of zeroth order and S^(t) is of first order in Bac- Thus the linearised 
Bloch equations become, for the zero-order and first-order contributions, respectively, to the S-components: 


dSl 

dt 

dt 

d^ 

dt 

d^ 

dt 


+ ^[S°At)-S°{oo)] =0 

- ujoSl + ^Si=0 
y T 2 

+ ujqS^ -I- 7 ^Sy — SaS° (t) coscot = 0 
-b 7^5'^ — SaS'^ {t) cos ujt — 0 


(25) 


where we introduced the resonance frequency ujq = jBz^dc = ~E/h, which depends on the level split¬ 
ting E, S = — poFcjCOS 9 which is a first-order term in the ac field 9 being the angle between 
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the ac field and the dipole moment. 5'°(0) is the initial value of 5'°(t) shortly after the field is applied, 
5 ^( 00 ) = —tanh (i?/2fcBr)/2 is the equilibrium value of the aligned spin, a = Ag/A, and we define 

X = H 1 — 4 ( 5 ° (00))^ /iksT. If one introduces raising and lowering operators = Si ± iS^ , then 
the equations for S'^ and S~ separate. The equations for S'^ becomes: 

dS+ (t) 


dt 


+ i ( wg — — ) 5'+ (t) — iaSSl (t) cos cot = 0 


(26) 


and the equations for S is the complex conjugate of the above. The solutions of these equations are given 
by the following expressions |56) : 

5 ° it) = S° (00) + [50 (0) - (00)] 

sx 

5'z (t) = - ^2 2 [coswt + Tiwsinwt] 

6 a [(wg — i/T2) cosujt — iwsinwt] 5 ° (00) 


5+ it) = 


(27) 


+ 


iLOQ-i/T 2 ) 

Sa [(wg + i/Ti — i/T 2 ) cosLot — iujsinujt] [S'° ( 0 ) — (00)] 

(Wg + i/Ti - i/T2)^ - 


the result for S~ it) being the complex conjugate of the above equation for S^it). The Bloch spins should 
now be related to the 2LS polarization in the electric field direction. The component py of the dipole 
moment along the direction of the electric field, in the diagonal basis as in Eq. (22), is now given by 
P|| = —(A(jj, + ^ax)pocos9 and using the average values of (ct^) and {ax) from the solutions of the Bloch 
equations one can obtain the dipole moment in the energy representation |56j : 


P\\ = -po cos 9 


f 2ASlit) AoiS+ + S-) 
\E E 


(28) 


Then, one must insert the deduced pseudo-spin values S^it), 5'+(t) and S it), Eqs. (27) to Eq. (28). Eq. (28) 
depends on the electric field linearly and can be easily differentiated with respect to the electric field, and 


this gives a formula for the complex dielectric constant e = 


dp II 
dF,., 


. For convenience, one may separate the 


F ^=0 


resulting formulae writing e = ic'jiEs6~^'REL)6~'i’i^'REs6~^'REL)- These are the real (e') and imaginary (e") parts 
of the dielectric constant. The imaginary part is interpreted as a dielectric loss (loss tangent, tan <5 = e"/e') - 
a parameter of the dielectric material that quantifies its inherent dissipation of electromagnetic energy (much 
like in a RLC circuit). One gets: 


e RES = 


Po‘^cos‘^9 

n 


Aq 

E 


(wg -|- w) F 


(^0 - 




tanh 


E 


(wg -|- Uj) T] 


, 1 + T ^2 i^O + <^) 
pg^COS^d f Ar 


(a;g - g;) T ^2 
1 -I- ~ w) 


e RES = 


n 


<-12 


E 


T. 


25° (0) -I- tanh 


To 


(—) 

\ 2 kBTi 


2 kBT 

-t/Ti 


1 -f (wg — U)) T 2 1 -|- (wg -|- Oj) T^ 
Ti 2 


tanh 


E \ 


1 -I- T^2 i^o + <kO) 


, pg^COS^0 

e REL = 


e REL = 


e REL — 



\ 2 kBT 


2 kBT, 

-tjTi 


2 kBT ) 

' 1 -f T2 w2 


I 1 

2 kBT) 

' 1 -f T^uj^ 

' E ' 

\ Tioj 


\ 2 kBT) l + T^uj^ 


(29) 


( 30 ) 
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where Tr 2 — T 2 In the adiabatic limit the initial value of the pseudo-spin is S'°(0) = — ^ ta nh(iii/2fcgT), 

when the 2LS energy eigenvalues are That makes the time dependent terms of Eq. (29) equal to 0 


shortly after applying the field. For an ensemble of 2LS, from the manipulation of the Bloch equations for the 
motion of the spatial components of a pseudo-spin 1/2 under periodic electric (and/or elastic) perturbations 
and taking into account the phonon relaxation mechanism, one then finds the explicit form of the expression 
for the dielectric constant m- 


e = iRES + srel 


1 


1 -I- iujT 


= (RES+iREL 


1 


1 -I- 


If^REL 


1 -I- 


(31) 


where we have put t = Ti. The above equation is the fundamental result of the low-frequency, relaxation¬ 
time approach to the complex dielectric constant of glasses. The typical energy splittings of the 2LSs in low 
temperature experiments correspond to frequencies in the range of « 10® Hz (GHz range, see Section 6), 
when the electric held frequency w is about 10® Hz. This justihes a low-frequency approximation a; <C wq. To 
obtain the resonant part we can also set — 0, which simplihes expressions (Eq. 29), remembering that 
E = Hujo- From the averaging over the dipole orientation angle 9 we get a prefactor 1/3: PqCOS^0 = ^Pq, 
where Pq is the conhgurationally averaged square 2LS electric-dipole moment. The real part of the relative 
dielectric constant for 2LS shows the temperature-dependent contributions (e(T) = e'(0) -I- Ae'(T), with 

|Ae'|«e'): 

A.' 


Ae' 


2LS RES 
1 


9 _ /\2 

= gP^^^tanh 


E 


2LS REL 


_)).^eosh 


2 kBT 
- 2 f E \ 


1 


\ 2 kBT J l+UJ^T^ 


(33) 


Eq. (32) corresponds to the resonant tunneling contribution to the dielectric constant, and Eq. (33) is the 


relaxational contribution. We neglect for now, for low w, the frequency dependence in the RES part so 
long as a; <C wq. The dielectric loss is described by the following formula (the resonant contribution being 
vanishingly small): 


A tan (5 


Ae" 


2LS REL 


2LS REL 


_ 


E \ 1 

2 kBT) 1+u^T^ 


(34) 


Integrating Eqs. (32 )-(34) over the parameter distribution of the 2LS and over the dipole orientation angle 0, 
using the expressions for the 2LS (E, r)-parameter distribution (see e.g. [M]), one can find the temperature- 
dependent contributions to the dielectric constant and dielectric loss: 


Ae' 

Ae' 


2RES 


2 Ppi r" 

3eoer y Ao,, 

TV, 


dE , f E \ 

— tanh - 

E \2kBT) 


' 1 - 


A, 


0,mm 

E 


i'E„ 


AtanJj 


2REL ^^O^rksT J 

pVo V 


dE I —^/l- 




.(E) ,_2 1 

r [2kBT) l+Lo^T^ 


(35) 


2REL 


SeoCrkRT Jao,, 


dE 


dE) 


'"V- 




E 


2kBT / 1 -h a;2T2 


If one extends, when appropriate and as a further approximation, the integration limits {Emax 00 and 
Ao.mm 0), then calculating (estimating, in fact) the E-integral one gets a characteristic logarithmic 
variation of the real part of the dielectric constant as a function of temperature: 


Ae' 


SeoEi- 


(*)• 


’'■Blinfj:' - A' 


2RES 


0 , 


T < 
T > 


2 kB 


(36) 


2 kE 


For ujTmin ^ 1 and at low temperatures, the contribution from relaxation to the real part as compared to 
the resonant contribution is negligible. Under the condition 1, however, the term l/r dominates in 

the r-integral and we obtain again a logarithmic variation with temperature: 


Ae' 

c' 


2REL 


V- 

IPPn 

SeoCr 


( To ) ’ 




(37) 
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A crossover between the resonant (low temperature) and relaxation (high T) regimes occurs at a characteristic 
temperature m 


To{uj) 


1 3 / OJ'Kph'^ 

+27?/u5 


(38) 


which for a thermal 2LS with E = A = ksT satisfies the condition ujt{Tq) = 1. The above equation arises 
from the celebrated expression (see e.g. [51 [57]) for the one-phonon relaxation time: 


T 


-1 



27: 


coth 



(39) 


arising from the longitudinal [1) and transverse (t) modes with deformation potentials 7 ;^* and sound speeds 
vi^ti respectively. In terms of {E,t) the 2LS distribution can be re-written as: 


ViE^T) 


P 

2r7r '^min {E)/t 


'^min (E) < T < 

^max (E) 


(40) 


(E) ~ E (sfefr) smallest (A = 0) and largest 


where Trmn{E) = ;^tanh(^27^) and 7max{^) = 

( 2 2 \ ~ 1 

^ -I- 2 ^ I is an elastic parameter of the solid 

and p is its density. 

Eqs. (36) and (37) show that with increasing temperature the T dependence changes from a decrease in the 
resonant regime to an increase in the relaxation one. At the temperature Tq{u}) there is a minimum. Thus, 
the sum of the two contributions has a characteristic V-shaped form, in a semi-logarithmic plot, with the 

/ y/3 

minimum occurring at a Tq roughly given by the condition ujTmin(kBTo) = 1, or ksTo^uj) = I ^ 70 ; ) . eoCr 


is here the bulk of the solid’s dielectric constant and we see that a - 2:1 characteristic behavior is predicted 
by the STM with the slope for T > Tq given by Eq. ( [^ . This behavior is indeed observed, but only in 
very pure (impurity-free) a-Si 02 [55]. However, in most multi-component glasses (chemically made up of 
GGFs as well as of GGFs, for example AlBaSiO (or BAS) glass: a-Al 2 03 -Ba 0 -Si 02 ), or for a contaminated 
mono-component glass, it is rather a V-shaped curve with a (roughly) -1:1 slope ratio that is often observed. 
This has been explained by our theory [19]. 


5.2 Results for the dielectric constant e' in a magnetic field 

We now derive the contribution to the dielectric constant e(uj) from the TWPs or ATSs sitting in the 
interstices between the RERs. One can treat the ATS again as an effective 2LS having lowest energy gap 
A£ = £i — £q = E = yjD'^ + Dqp'^ for “weak” fields. Within this picture, the linear-response quasi-static 
resonant and relaxational contributions to the polarizability tensor are extracted according to the general 
2LS approach described in the previous Section 5.1, to get [45115^ 

(41) 


and 


where 


a 


REL 


AkeT . 


dE 



E^E, 


PifiPji/ 


cosh 


E 


2 kBT 


5{E-A£) 


;Pi 


= E 


Pi^zPl, 




EiEj 

E^ 


PifiPj i 


(42) 


(43) 


contains the single-well dipoles Pi = qsn. This expression assumes vanishing electric fields and no TS-TS 
interactions, a situation which does not wholly apply to the experiments. To keep the theory simple one can 


still use Eq. (41) and the analogous one for the relaxational contribution to the polarizability. Eq. (41) must 


be averaged over the random energies’ distribution ( 10 ) ([.. .]a«, responsible for the high sensitivity to weak 
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fields) and over the dipoles’ orientations and strengths ((...)). For a collection of ATS with riu, > 2 wells 
this averaging presents serious difficulties and one must resort to the decoupling: 


g^,5{E - A5) ~ ■ 5{E - A£), 


(44) 


where [S{E — AS)]av = gATs{E, B) is the fully-averaged density of states. To calculate one can envisage 
a fully isotropic distribution of planar n^,-polygons to obtain [5^ : 


gtiv = 


- 1 


-^{n^ - 2)E'^ + Dlip' 

P^ ^2 




(45) 


The second term in the numerator of Eq. (451 gives rise to a peak in Se/e at very low B, while the 
first term (present only if > 2) gives rise to a negative contribution to 5e/e at larger B which can 
win over the enhancement term for all values of B if Domax ^ (Domin, Domax corresponding to 

cutoffs in the distribution of ATS energy barriers). The observations in Duran and BK7 indeed show a 
significant depression of e'{B) for weak helds [7], thus giving direct evidence for the existence of ATSs with 
riu, > 2 in the multisilicate glasses. Carrying out the averaging [... ]av one gets analytical expressions for 
the polarizability; the uniform average over orientation angles 9 must be performed numerically (although a 

i)- 


very good approximation is the replacement of with P(p'^ in the averaged expression, corresponding to the 
replacement cos^ 9 


The expression for the resonant part of the polarizability Eq. (41) should be averaged over the probability 
distribution of parameters Eq. ( |10[ ). Averaging over the probability distribution 

/ dEidE 2 dE^V'XTS -^o) d ^i) can be done as in the previous Sections, to get to the 

expressions 


_ poo -I / 7? \ 

»res{T. B) / <i£^ta.,h (—) 


o„ (T, B) =v>%P-pI I gtanh 
ai{T,B)= 1 p-^^ r'Etsmhf ® '' 


f """ [E- _ 


D JDo 


xS{E-AE)= ai (T, B) + ao {T, B) 


E 


\ 2 kBT) 


dD 


l•D„ 


— dDoDoS (E - A£) 

■OO ^Dorr^a.. 


^ ^min -^0 


(46) 




a 


0 


6 {E- A£) 


with the energy gap A£ = y/D'^ + 

Similarly to the calculation for the density of states Eq. (13), the integrals in Eq. (46) can be reduced as 
follows IMI: 


ao(r,S) 


. f E \ 

^ tanh - 

E2 \ 2 kBT) 


/“ f i».h j 


'Dorr 

pDor 


y/E^-Dl 


DodOg 


E2 - 
Dq dDo 




rEc2 


{2kBT) jBorrrrrr 

^tanh ^In ( 

E^^^^^[ 2 kBT) 2 ^\ ) 

/ \ 1 / 7:^2 7-42 , _2 \ H 


(47) 


l^Eci \ ^ / \ r 

Je.,E^^ ^[2kBT)2^ 
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Material and Temperature 

'KXATsP*plle-re^o 

D K 

-E^min ?Av 

Dormn i 

Domax i 5a,K2 

BK7 15 mK 

0.089 • 10"'’ 

0.03 

1.668 • 10" 

4.576 • 10'’ 

Duran 15 mK 

0.052 • 10"'’ 

0.021 

2.457 • lO'^ 

4.151 • lO'^ 

AlBaSiO 50 mK 

0.89 • 10-'> 

0.015 

2.440 • lO'^ 

3.080 • 10'’ 

AlBaSiO 94 mK 

3.75 • 10-'> 

0.025 

1.225 • lO'^ 

1.589 • 10'’ 

AlBaSiO 120 mK 

3.09 • lO"'^ 

0.0227 

1.767 • lO'^ 

2.248 • lO'^ 


Table 5: Fitting parameters for the dielectric constant in a magnetic field for three different types of glasses. 


Temperature 

1 TXATsP*p\/f^rf -0 

D K 

Dormn ? 

Doma. i K2 

50 mK 

4.38 • lO"'’ 

0.015 

0.076 • 10^ 

3.047-10^ 

70 mK 

12.22 • 10"^ 

0.0486 

0.600 • 10^ 

2.662 • 10^ 

100 mK 

13.63-10-^ 

0.0486 

3.035 • 10^ 

7.616 • 10^ 


Table 6: Fitting parameters for the Si02+a;CyH2 glass for different temperatures. 
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with Eci,Ec 2 as in Section 4.1. For the whole mass of the glass: 

- ^ (ai (T, B) + ao (T, B) - a {T, 0)) = Ao = xats^d 
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where xats is the volume concentration of ATS and a {T, 0) is: 
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(48) 
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(50) 


(51) 


The relative change of the dielectric constant is expressed by Eq. (47), Eq. (48) and Eq. (50): 

Ae'(T,B) Aa {T,B) xats , u\ , irr b,\ irr a\\ 

--= Xats -=-(ai (T, B) + og [T, B) - a (T, 0)) 

e CoSr eo^r 

Eq. (51) with Eqs. (47)-(48) describes well the experimental data for different glasses, as is shown in Figs. 8p0 
and with the fitting parameters presented in Tables Hi 

For the sake of clarity, the data and curves in Figs. and have been shifted apart vertically. 


5.3 Results for the dielectric loss in a magnetic field 

The dielectric loss (or loss angle 6 ) for a dielectric substance is a measure of the power lost in dissipation 
and is obtained as the following expression: 

e" e" 

tan 5 = — = - (52) 

e' eoCr- 
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Figure 8: (colour online) The relative dielectric constant variation as a function of the magnetic field for 
AlBaSiO (BAS) glass [7]b, BK7 [7] and Duran Hb glasses. With best-fit parameters as in Table the 
curves are the results of our theory in the “weak field” approximation with (and, for AlBaSiO, without) 
higher order correction. From |53j . 



Figure 9: (colour online) Relative dielectric constant variation as a function of the magnetic field and tem¬ 
perature for AlBaSiO (BAS) glass [7]b. With fitting parameters as in Table the curves are the result of 
our theory in the “weak held” approximation. From |53j . 
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Figure 10: (colour online) Dielectric constant variation as a function of the magnetic field and temperature 
for the Si02+a;CyH2 glass [50]. Fitting parameters as in TableFrom [55] . 


where e" = enEL is the imaginary part of the dielectric constant, typically evaluated in the relaxation 

time approximation (Section 5.1). It should be pointed out that the 2LS STM does not describe well the 
temperature- and frequency-dependence of e" in glasses. The reason is that the theory works well only in 
the low-frequency regime. Since e" (w) should be linked to the real part e' of the dielectric constant through 
the Kramers-Kroning relation, one can see that the knowledge of the low-frequency behaviour of e' {uj) is not 
enough to reproduce the correct form of e" {uj,T). However, inclusion of the ATS contribution does seem to 
improve the agreement between theory and experiment [53] , as least in the case of BK7 (optical glass) which 
is probably the best approximation to a fully networked glass. 

The relaxation part of the dielectric constant depends on the volume concentration xats = f^ATSP of ATS’s 
and is expressed by: 

^REL — 1 = Xats chhel (53) 


where a eel should be averaged over all parameters by means of the distribution function Eq. (lOl; following 
similar calculations as in the previous Sections we get: 
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(54) 


The relaxation time for ATSs at low temperature and in a magnetic held is now found (after a long calculation 
that will be reported elsewhere) to be given by the following expression |5^ : 


‘ATS 


= T [E,if) = 


E^ 


ID^) 


E^ {E^ - iD'^) 

rtanh(2fcfT) rtanh(2l|T) 


= T 


-1 


(A,D) 


(55) 


where as usual the A-B phase ip is directly proportional to the magnetic held B. It appears, therefore, that 
the total dielectric relaxation time, obtained through its inverse: 


1 _ 1 1 
Ttot T2LS TATsi^f) 


(56) 


must diminish in a non-trivial manner as the magnetic held is switched on. This very interesting prediction 
of the present theory appears to be conhrmed explicitly, albeit only qualitatively, in the laboratory and for 
some multi-silicate glasses so far only via the work of a Russian group at liquid-He temperatures EH. A 
systematic study of the magnetic-held dependence of Ttot in the multi-component glasses is still lacking. 


20 


















Temperature 

TTXATsP*Pi/e-PO 

Dmin 7 K 

Dormu i 5a, K2 

Domax i 5a, K2 

r=7'fc|, (sK")"i 

AlBaSiO 

77 mK 

1.5410"'’ 

0.0209 

1.9810'' 

4.9610" 

5.010'' 

88 mK 

1.3510"'’ 

0.0206 

1.9810" 

4.9610" 

4.010" 

96 mK 

1.1010"'’ 

0.0213 

1.9810" 

4.9610" 

4.410" 

BK7 

15 mK 

2.0210"'’ 

0.0287 

0.6910" 

0.6610" 

3.3410" 


Table 7: Fitting parameters for the dielectric loss in a magnetic field in the AlBaSiO (BAS) and BK7 
glasses. 


The probability distribution function for the ATS dielectric relaxation times turns out to be rather different 
from that of the standard 2LS case, one finds indeed 1551: 


P(A,r) = 




bP (t 'Train) ij'max 't) 


(57) 


(with suitable (^-dependent boundaries) where (with T a new ATS-related elastic constant). Tmax{E) = 


/ r tanh (^ 2 fchT j maximum ATS relaxation time and the minimum allowed is Tmin{E) = (5/6) Traax{E). 

The ATS relaxation-time distribution is therefore very narrow-ranged in r and also very singular (albeit in 
the B = 0 case only). This has important experimental consequences that will be discussed elsewhere. 

It is however not convenient to switch to r as an integration variable. Then one writes, integrating over Dq 
first via: 
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(58) 


Thus the integral in Eq. (541 becomes: 
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The integral in Eq. (59) has two special points Di (E) = ^ E'^ — and D 2 (E) = ^E'^ — 

hence Di {E) < D 2 JE). The 0-conditions divide the integral in Eq. ( [5^ into two terms with different 

intervals in the energy value: E^i < E < Ec 2 and E > Ec 2 , and the integal for E < E^i vanishes: 
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(60) 


dielectric loss variation in a magnetic field: 


Substituting Eq. ( |55[ ), Eq. ( |60[ ) into Eq. ( |53[ ) and then into Eq. (52) we obtain the final formula for the 
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(61) 


withe Ed and E^ as in previous Sections. The fitting of relative dielectric loss variation in a magnetic 
field is shown in Figure El using the best-fit parameters from Table One can see that, once again, the 
experimental data are very well reproduced by the present theory and with fitting parameters very similar 
to those extracted from the study of the (real part of the) dielectric constant (Section 5.2). 
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Figure 11: (colour online) The relative dielectric loss as a function of the magnetic field and temperature in 
the AlBaSiO (BAS) (a) and BK7 (b) glasses (data from [7]b). The continuous curves are from the present 
theory. In the inset of (a) we show that a faint peak seen experimentally at very weak fields can also be 
explained by the theory. From [53) . 
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6 The Magnetic Field Dependent Polarization Echo Amplitude 


6.1 The Polarization Echo Experiment 

The experimental detection of electric and phonon echoes in glasses is one strong convincing argument for the 
2LSs’ existence. Echoes in glasses are similar to other echo phenomena such as spin echo, photon echo and so 
on. But only at very low temperatures the relaxation of the TSs becomes so slow that coherent phenomena 
like polarization echoes become observable in the insulating glasses. 

The essence of the effect is the following (see Fig. 12). A glass sample placed in a reentrant resonating 


cavity (“Topfkreisresonator”) is subjected to two short ac electromagnetic pulses at the frequency of about 
1 GHz separated by a time interval ti 2 . The duration ti and T 2 of these pulses should be much shorter that 
all relaxation processes in the observed system. The macroscopic polarization produced by the first pulse 
vanishes rapidly due to the distribution of parameters of the TSs in glasses. This phenomenon is similar 
to the well-known free-induction decay observed in nuclear magnetic resonance (NMR) experiments. The 
“phase” (energy-level populations) of each TS develops freely between the two exciting pulses. The second 
pulse causes an effective time reversal for the development of the phase of the TSs. The initial macroscopic 
polarization of the glass is recovered roughly at a time ri 2 after the second pulse. Since the thermal relaxation 
processes and (see later) spectral diffusion are strongly temperature dependent, polarization echoes in glasses 
can be observed in practice only at very low temperatures, typically below 100 mK. The echo amplitude is 
proportional to the number of TSs that are in or near resonance with the exciting microwave pulse and that 
do not loose their phase coherence during the time 2 ti2 [H]. It should be pointed out that, due to the wide 


I pulse II pulse Echo 


1 



a 

, 






r T,2 ' 

Ti2 1 




Polarization echo experiment 



(a) 


(b) 


Figure 12: (colour online) The two-pulse polarization echo experiment. Hahn’s vector interpretation on the 
right hand side is for NMR’s spin-echo experiment. 

distribution for the parameters of the two-level systems in glasses, the description of polarization echoes in 
glasses is much more complicated than in the case of nuclear spin systems. In analogy to the two-pulse echo 
in magnetic resonance experiments this phenomenon is referred to as the spontaneous echo. 

The polarization echo phenomenon can help to understand more about the microscopic structure of TSs 
in glasses and gives different kinds of information. The analysis of these experiments follows that for the 
equivalent magnetic case, except that the TS problem is complicated by three factors. First, the elastic or 
electric dipoles are not aligned with respect to the driving field and a calculation of the echo signal involves 
an average over their orientations. Secondly, for a given pumping frequency lo there exists a distribution of 
induced moments (electric or elastic) and relaxation times, which should be included in the analysis. Finally, 
in electric echo experiments the local field seen by the TSs is not equal to the applied field, and a local-field 
correction factor must be used when evaluating absolute values of the dipole moment [42] . 

In the polarization echo experiments at radio frequencies and at very low temperature of about 10 to 100 mK 
it has been shown that the TSs in glasses couple directly to the magnetic field [SU |S2|. Unexpectedly, the 
amplitude of two-pulse echoes in the BAS glass was found to be strongly dependent on the applied magnetic 
field showing a non-monotonic (even oscillatory) field variation. In subsequent papers [HI , such behavior 
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was attributed to the existence of nuclear electric quadrupole moments (NEQM) for some tunneling particles 
(having nuclear spin / > ^) interacting with the magnetic field and with gradients of the internal microscopic 
electric field. The NEQM model is based on the consideration that the levels of tunneling particles with non¬ 
zero nuclear quadrupole moment exhibit a quadrupole splitting, which is different in the ground state and in 
the excited state of a tunneling 2LS. The magnetic field causes an additional Zeeman splitting of these levels 
giving rise to interference effects. In turn, these effects cause the non-monotonic magnetic field variation of 
the echo amplitude. 

The amplitude (or integrated amplitude) of two-pulse polarization echoes of four types of silicate glasses is 
shown in Fig. |13[ a) as a function of magnetic field |62j . In contrast to many other low-temperature properties 
of glasses the influence of the magnetic field on the amplitude of spontaneous echoes is obviously not universal. 
BK7 and Duran show similar effects, although the concentration of magnetic impurities differs by at least a 
factor of 20. Perhaps a most remarkable result of the measurements is the fact that Suprasil I (very pure 
a-Si02) shows no measurable magnetic field effect. While Duran, BAS and BK7 contain nuclei with non-zero 
nuclear quadrupole moment, Suprasi I is virtually free of such nuclei. This fact is used to provide justihcation 
for the nuclear quadrupole model. The variation of the echo amplitude with the applied magnetic field is 
similar for Duran, BK7 and BAS, but not identical. All three samples exhibit a principal maximum at very 
weak fields, B 10 mT, but only BK7 has a relevant second maximum and a hint to an oscillation in B. At 
high fields the amplitude of the echo rises well above its value at zero magnetic field and seemingly saturates 


(yet, see Fig. 


In Figure 13 


this is very similar to what happens to the inverted heat capacity, —Cp, as a function of B). 
the amplitude of spontaneous echoes in the BK7 glass is shown as a function of the applied 



Figure 13: a) The integrated echo amplitude as a function of the magnetic held for different silicate glasses: 
BK7, Duran, AlBaSiO (BAS) and Suprasil I. All data were taken at T=12 mK, ti2=2 ms, and roughly 1 
GHz, except for Duran, where the delay time was ti2 =I.7 ms [52]- b) The amplitude of two-pulse echoes in 
BK7 glass as a function of the magnetic Held for different values of the waiting time ti 2 between pulses. All 
data sets were taken at 4.6 GHz and 12 mK except that for ri2=2 ms which was taken at 0.9 GHz. 

magnetic Held for different delay times ti 2 between the exciting pulses. We can see obvious differences for 
different values of ti 2 and that a second maximum (the “oscillation”) is not always present. These Hndings 
necessitate a good theory for spectral diffusion in real glasses and this theory remains to be accomplished. 
The most remarkable fact about these experiments on echoes from glasses in a weak magnetic Held is that 
the strong magnetic effect is not conHned to the inorganic, silicate glasses. Figure shows the amplitude 
of spontaneous echoes in partially deuterated and in ordinary amorphous glycerol as a function of the weak 
magnetic Held B. In the case of ordinary glycerol (do) there is very small change of the echo amplitude with 
B. However, for partially deuterated glycerol (da) a change is much more noticeable, of a different shape and 
duration. This experiment seemingly provides proof that the magnetic effect is of nuclear origin, for the two 
amorphous glycerol samples differ in the content of nuclei carrying a NEQM. Glycerol-dg has none, other 
than the natural abundance of deuterium, some 125 ppm, and of about 500 ppm, concentrations which 
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(a) (b) 

Figure 14: (colour online) The integrated echo amplitude as a function of the magnetic field at r=13 mK, 
generated in partially deuterated glycerol (glycerol-da) and in ordinary glycerol (glycerol-dg). On the left- 
hand side the figure shows that the echo amplitude for deuterated glycerol-da is much more sensitive to the 
magnetic field in comparison with non-deuterated glycerol-do, shown at the right-hand side (from [5]). 


are however a factor 10 too weak to account for the observed magnetic effect. Glycerol-da contains instead 
37.5% D (1=1) and 62.5% H (7=1/2). However, glycerol-dg (nominally 100% D) displays a magnetic effect 
that is only 10% larger than in glycerol-da while glycerol-ds (62.5% D) displays a smaller magnetic effect 
than glycerol-da [HISS]- All this hints to the fact that the effect does not scale with NEQM concentration. 
As reproduced in Fig. the integrated echo amplitude as a function of the waiting time ri 2 in amorphous 
partially deuterated glycerol-ds (that is CaOaHaDs instead of ordinary CaOaHg) shows exponential decay 
with spectacular oscillations at zero applied magnetic field. On the other hand, all oscillations disappear for 
the relatively weak magnetic field of 150 mT [SIIMIES]. These findings are extraordinary, especially when 
combined with the observation that for all multi-silicate glasses the oscillations of the echo amplitude in ri 2 
are absent for all values of the magnetic field (a fact that the NEQM approach cannot explain). All these 



Figure 15: (colour online) The integrated echo amplitude as a function of the waiting time ti 2 at zero 
magnetic field (red curve) and in weak magnetic field of 150 mT as generated in deuterated glycerol-ds (from 

[MIISS])- 


findings are collectively hard to explain on the basis of the 2LS STM or starting from more microscopic 
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models and so far only the NEQM has been able to provide the beginning of a rationale for some of these 
startling experimental results. 


6.2 The polarization echo in a magnetic field: Schrodinger eqnation formalism 


Using the density matrix formalism of quantum mechanics, it is possible to obtain, for a collection of 2LSs 
only, the following expression for the average dipole moment near t = 2 ti2'. |53j 
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(62) 
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(64) 


where: F = Fo(t) cos(a;t) is the pulsed external electric field (having intensity Fq and MW frequency w), 
Po is the local bare electric 2LS dipole moment, y = cos0, 0 is the dipole’s orientation angle wrt Fq, 
Q,fi = AqPo • Fq/ HE is the Rabi frequency of the resonant 2LS, flc = + (w — E/hY is the generalized 

(non-resonant) Rabi frequency, z = E/h, — to and where 7 = 2 t~^ is the 2LS intrinsic phonon damping (t 
being the 2LS phonon-induced relaxation time (Section 5.1)). We have also defined a fluctuating 2LS energy 
gap E(t) = E + tiAuj{t) with hAuj{t) the contribution due to interactions with the surrounding thermal 2LSs, 
and a service function: 


s{t) 


-1-1 if 0 < t < ri2 
-1 if t > Ti 2 


(65) 


This result reproduces (and improves) the derivation by Gurevich et al. [B5]. It should be stressed that the 
(normalized in ojq = E/h) spectral function a {ojq) = is in practice, for <C w, very narrowly 

centered at the pumping frequency w and can therefore be replaced by a Dirac’s S (wq — w), thus resulting 
in the strictly-resonant approximation which is normally employed. Things change considerably, however, 
when the density of states g{E) is non-uniform, as is the case for the ATS model. 

The above result for the echo signal from a collection of 2LS can also be obtained - and from first principles 
- from a lengthy but straightforward Schrodinger equation treatment in which high-frequency modes are 
neglected and phonon-damping is treated in a phenomenological way. In the most rigorous way, one obtains 
for the echo signal: 

p|| (t) = p+ (t) cos (wAr) — p_ (t) sin (wAr) (66) 

where At = T2 — ti (whence wAr typically as large as 10^) and where: 
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(67) 


in which 2 ti2 = ti + T12 -f T2 -l- T12 ~ 2 ti 2 is the total elapsed time at the echo signal’s centre. The above 
Eq. (66) gives the expectation value of the polarization per TS in the direction of the applied electric field. It 
still needs to be averaged wrt all STM parameter distribution and over a uniform orientational distribution 
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of 2LS dipoles po- \=t is again the phonon relaxation rate. If wAr is neglected, then Eq. (64) (after 


averaging) is recovered for Py = pj[. 

We are now in the position to extend the polarization echo’s calculation to the case of the ATS model 
describing glasses in a magnetic field; the point of view will be taken that a background of ordinary 2LS’s - 
insensitive to the magnetic field - also exists in the glass. This is in line with our generic cellular model for 
the real glasses (Section 2). 

One starts with a collection of 3LS {n^ = 3 is not only computationally convenient, but physically correct 
as explained in Section 2), but with the single ATS Hamiltonian written in the energy representation: 
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( 68 ) 


where the diagonalizing matrix S = Si^p) is magnetic-field dependent, the Ei are the (H-dependent) ATS 
energy levels and the Pi are the wells’ electric dipoles. As in the treatment of Gurevich tt al. [56] there is also 
a phonon bath, but this will be treated - as always - phenomenologically and resulting in a phonon-damping 


exponential. The second term in Eq. (68) causes irrelevant energy-level shifts and produces an extra matrix 


term AH'(t) = (Ay) of which the only relevant element (see below) is 


Aoi — A 


10 


= ^ -Pfc • Fo Sok(<p)S);k(:^) cosut 


(69) 


fc=i 


The Ay cause transitions between the ATS levels 0, 1, 2 when the pulses are applied. In the weak magnetic 
field limit (most appropriate for the echo experiments) and in the approximation D <C Dq that we always 
use (and that is always confirmed by out best fits to the data), one quickly discovers that the second excited 
level remains unperturbed and one can make use of the “effective 2LS approximation” (where, however, 
the ground-state wavefunctions of the three wells mix). One can then repeat the Schrodinger equation (or 
density-matrix, for that matter) calculation carried out for the 2LS case, at the cost of introducing a complex 
Rabi frequency: 

^0 = ^ (70) 

n 

The evolution of the generic ATS during and in the absence of pulses can then be followed exactly, in much 
the same way as before, except that in order to simplify the formalism it is convenient to introduce from the 
outset an orientationally-averaged Rabi frequency (now a real quantity): 

(71) 

the bar denoting the average wrt 3LS base-triangle’s orientations. Replacing Oq with Eln before carrying out 
the averaging of the sample’s polarization is our main approximation, allowing for a considerably simplified 
treatment and leading to the magnetic-field dependent expression [55] : 
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PiFo ^ 


h 




(72) 


Here, pi is a single-well (averaged) electric dipole and E = hujo = is the usual magnetic- 

field dependent lower energy gap in the weak field approximation. The above approximation for Hq treats 
incorrectly the ATS’s that have Fg roughly orthogonal to the ATS base triangle; luckily these have Hg « 0 
and do not contribute to the echo signal. 


Proceeding as for the derivation of Eq. (661 one finds that there is a magnetic contribution to the (partly 
averaged) polarization of the sample from the generic ATS given by [59] : 
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Now, \ = T is the magnetic ATS phonon relaxation rate given by Eq. (54), the generalized Rabi frequency 


is again given by Vlc=\J n|j+(a;o—and: 

$ {t) = ujoit — 2ri2) + ujAt 


(74) 


is the appropriate time argument. From this, it is obvious that the time at which all ATS (regardless of their 
energy gap E = huJo) will be refocused is t = 2 ti2 and this determines the echo’s peak position (if the echo 
signal has a reasonable shape, which is not always true |52j'l. The measured echo amplitude’s contribution 
from the magnetic ATS is therefore (allowing for an arbitrary amplification factor Aq): 
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where d is the sample’s thickness, 0(T>, Do) is the usual theta-function restriction for the integration domain 
(previous Sections) and where a final orientational averaging wrt the angle j3 = BSa (defining the A-B 
phase p, see Eq. ([^) is in order. At this point one deals with the delta-function’s constraint and the energy 
parameters integrations in the usual way, to arrive at, after a lengthy calculation |59j : 
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where we have defined the functions: 
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with 01^2 = ^gti ,2 the so-called pulse areas. 2 are as in the previous Sections, whilst D 12 (p) = 




- Pomax,m^nV’‘^ and E = fu^o- 


In going from Eq. (73) to Eq. 


(|76[) we have tacitly made some assumption on the (fully averaged) spectral 


diffusion term Au){t)s{tyit ^ theory of spectral diffusion (SD) for the magnetic multi-welled 

ATS is a chapter still open, however we can safely assume that what was found by many Authors for NMR’s 
spin-echoes and for the 2LS polarization echoes in glasses holds for the ATS as well. Namely, that there is 
a wide range of waiting times where the decay of the echo amplitude is a simple exponential in ri 2 so that 
one can replace the SD term with ^ where TipiT) is a SD characteristic time depending only on 

temperature. There must be a SD time t^( 3 ) for the ATSs as well as a SD time r,^( 2 ) for the standard 2LSs’ 
ensemble. For the latter, theory shows j67l [68] that this parameter is independent of the energy gap E and 
thus for the ATS we shall assume the same and, moreover, that (like for the phonon damping rate and Rabi 
frequency) its dependence on the magnetic field is weak or absent. This allows us to lump the SD problem 
together with phonon damping, yielding an overall exponential relaxation rate: 


w{E,D) = t^-^ +t-\E,D) 


(78) 


in which the SD time is typically much shorter than the phonon-damping time t and depends only on 
temperature through: 


= catsT 


(79) 
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BAS (a-Ba0-Al20j-Si02): Echo Amplitude vs B 


T=12.0 mK (sample 1), 20.0 mK & 56.0 mK (sample 2) 



Figure 16: (colour online) Magnetic field dependence of the polarization echo amplitude (relative to its 
value at “high” fields where saturation occurs) for the AlBaSiO glass [22] (also referred to as BAS) at 
given experimental conditions. We believe two separate samples have been used. Continuous curves from 
our theory. Nominal frequency 1 GHz, T2 = 2 ti— 0.2 us. Inset: behaviour of the ATS DOS for the same 
parameters (the physical origin of the effect). 


with Cats an appropriate constant. The assumption of an overall simple-exponential decay of the echo 
amplitude with ti 2 is well verified experimentally m- 

We now make use of Eq. (761 to fit the experimental data for the multi-silicates, the idea being that the 
total amplitude is A {(p) = A2LS + AA (tp) (which must be averaged wrt the ATS magnetic orientation angle 
/3). Fig. 16 shows the experimental results for the relative echo amplitude in AlBaSiO (BAS glass) as a 
function of the magnetic field; values of B up to 0.6 T have been explored and for three temperatures. The 
data are fitted with our theory with parameters as reported in Table [^ The agreement between theory and 
experiment is highly satisfactory, given the simplifications used in the theory. There is only one minimum 
in A{B) and the inset in Fig. shows that again it is the ATS density of states (DOS) that is responsible 
for the magnetic effect (Section 3). Indeed, by enforcing the strict-resonance condition a{E) —>■ S{E — fuE) 
Eq. (761 collapses to a quantity very much like the DOS (convoluted with slow-varying corrections) and 
with the same behaviour, thus reproducing the main shape of AA(i?). It is the non-resonant convolution of 
this quasi-DOS with other E-dependent functions that produces the rounding of the minimum and the B~^ 
saturation. Interestingly, though t, the phonon-damping term plays a main role in the rounding of the 
high-E tail to a B~^ (as observed) saturation. The ATS approach predicts also a linear in B intermediate 
decay regime of the echo amplitude, and this is often experimentally observed. Next, in Fig. 17 we present 
the comparison of theory and experiment for data for the echo amplitude in BK7 (good optical glass, hence 
devoid of true microcrystals, bust nevertheless containing RERs) at two different values of the waiting time 
ti 2 . It is remarkable how our theory, despite the simplifications and the total absence of multi-level physics 
(as advocated by the NEQM approach), can reproduce all the features of the experimenatl data, including 
every change of curvature in A(B) vs. B. A rough fit, not aiming at high agreement, reproduces the 
two maxima (and minima) that the NEQM approach takes as indication of the multiple (rapid) oscillations 
ensuing from the quantum beatings due to the Zeeman- and NEQM-splitting of the generic 2LS |23| ■ There 
are never more than two observed minima, in practice, and these can be reproduced by our simple ATS 
model. Finally, in the inset of Fig. |17[ we show what the experimentalists missed by not exploring higher 
magnetic-field values. Using the simple-minded correction for the lower energy gap at higher fields, we plot 
the expected behaviour of A(B) for intermediate fields. After the two minima, there is only an apparent 
saturation and new interesting features should characterise A(B) at higher fields (B >600 mT), just like it 
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BK7: Echo Amplitude vs B (T^12 mK) 



Figure 17: (colour online) Magnetic field dependence of the polarization echo amplitude for the BK7 glass 
[62] at given experimental conditions. Dashed curves (rough fit) and continuous curves from our theory; 
there are no more than two observable maxima or minima (no true oscillations). Nominal frequency 0.9 
GHz, T 2 = 2ri=0.2 /xs. Inset: our prediction for the higher magnetic field regime {B* as defined in Section 
5). 


Glass type 

Drain 

(mK) 

Domin 1 g 1 

(kA^) 

Domax 1 g 1 ‘5'a 
(KA 2 ) 

F-i 

{usK^y^ 

CATS — C2LS 

PiFo 

D kV m-i 

tan w At 

AlBaSiO 
(sample 1 ) 

17.74 

0.95x10^ 

2.13x10^ 

9.22x10® 

5.008 

0.461 

0.247 

AlBaSiO 
(sample 2 ) 

27.20 

1.14x10^ 

8.96x10^ 

2.57x10® 

3.825 

0.450 

0.245 

BK7 
(1.5 ^s) 

16.76 

0.92x10^ 

1.34x10^ 

8.91x10® 

1.03 (*) 

0.60 

0.207 

BK7 

(6 ^is) 

15.94 

0.89x10^ 

3.31x10^ 

3.25x10® 

5.72 (*) 

0.98 

0.204 


Table 8 : Fitting parameters for the echo amplitude’s magnetic field dependence. (*) For BK7 (best-fit 
parameters only), cats only is involved. 


happens for the dielectric constant (Section 5). A full description of the effect, however, requires a calculation 
involving all three ATS energy levels. 

6.3 Amorphous glycerol and the so-called isotope effect 

We now come to the astonishing case of amorphous glycerol (C 3 O 3 H 8 vitrifies around Tg ~47 K), deuterated 
and natural. The first question is whether our model applies to this system, which nominally is single¬ 
component. We firmly believe it does and therefore that the a-glycerol polarization echo experiments can 
also be explained by the presence of RERs/microphasing in the samples. 

In Brandt’s Ph.D. dissertation it is reported that the liquid glycerol, from which the glass samples 
were made of, were contaminated by water (see Table |^. Moreover, experiments on samples with different 
deuterium molar content and different before-cooling open-air shelf-storage times gave definitely different 
results. The available experimental data are reported in Fig. |18[ One observes a much greater variation in 
the experimental data for the C 3 O 3 D 3 H 5 samples left in the air before freezing than in the similiarly prepared 
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Sample 

Ghemical purity 

Water content 

glycerol-do 

99.9% 

ppm H 2 O 

glycerol-ds 

99% 

= 1.5% H 2 O 

glycerol-ds 

98% 

= 0 . 11 % H 2 O 

glycerol-ds 

98% 

no information 


Table 9: Purity and water-contents data for the studied glycerol samples. 




Pulsabstand (fis) 


(a) 


(b) 


Figure 18: (colour online) Relative integrated echo amplitude vs. T 12 for various deuterated amorphous 
glycerol samples at R = 0, nominal frequency 0.85 GHz, T =13.5 mK. The samples termed “old” (alt) were 
left liquid several days in the air before cooling; two separate fresh preparations (neu(l) and neu( 2 )) were 
also employed, obtained from producer-sealed containers just before cooling. From [69]. 


C 3 O 3 D 5 H 3 samples. 

Glycerol is a highly hygroscopic substance, owing to the polarity (similar to H 20 ’s) of its molecules which tend 
to shed one of the three hydrogens (or deuteria) attached to the oxygens, but none of those firmly attached 
to the carbon atoms. The equilibrium process is of the type, e.g., C 3 O 3 D 3 H 5 <—^ G 303 D 2 H^-|-D+ causing 
similar equilibrium concentrations as in the familiar dissociation process H 2 O —>■ OH~-|-H+. Thus it is not 
unreasonable to imagine that the water contained in the air causes some hydrogens to substitute for some 
0-attached deuteria, in time, when absorbed in liquid glycerol-d 3 , - but not in the case of glycerol-da. The 
glycerol-da hygroscopic effect on the echo would support the NEQM theory, but the (albeit smaller) effect in 
glycerol-ds does not. Moreover, oddly enough the more hygroscopically contaminated samples (“alt”) yield 
the smaller effect on the echo’s amplitude oscillations, whilst the fresher samples (“neu” ( 1 ) and ( 2 ) give 
wider and different (neu(l) 7 ^ neu( 2 ) amplitude oscillations. 

The lack of reproducibility of the effect leads us to believe that the origin of the observed magnetic phenomena 
lies again in the formation of RERs/microcrystallites during the freezing process, which in glycerol are 
nucleated by the presence of dissociated H 2 O. The concentration of micro- or nano-crystals is therefore very 
much sample-dependent. 

We remark from further tests [53] that (Section 9) the number of coherently tunneling ions (most probably 
D'b in glycerol-ds) may be large and/or that the electric-field intensity (unknown for the glycerol experiments) 
may also be huge (thus renormalizing the value of Dmin way upwards). 

The lack of reliable experimental dielectric data for amorphous glycerol makes the htting of the polarization- 
echo data frustrating. All the more so, when considering that - so far - no evidence for a magnetocapacitance 
effect in a-glycerol has been found (albeit at the nominal frequency of 0.986 GHz) [69]. This last fact, however, 
could be explained by the suspected large value of Dmin, for the order of magnitude of the magnetocapacitance 
is given, in our theory, by the combination t:P*p\xats /(Section 5) and for a-glycerol Sj. = 42.5 
and Drain are much larger than for the multi-silicates. Low-frequency (kHz range) and weak electric-field 
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intensity measurements of the dielectric constant of a-glycerol would be most useful. These lacking, we 
must limit our discussion of the magnetic dipole echo phenomena in a-glycerol to the qualitative level of 
understanding. 

The integrated echo amplitude, lEA, is readily obtained from Eq. (73), after the averaging procedure as 
in Eq. (75) and a subtraction of the pump-frequency mode cu are carried out (rotating reference frame 
description): 


A TT ^ \ A d, At:H^P 

AIEA{(p) =-Ao - XATS — 

SoSr ro 

dD , 
—— tanh 


E 


I D„ 


D 


E 


E^ 


2 kBT I E‘^-D^ 


nla{E) 


2w 


+ {UJQ — oj) 


(80) 


X [S ( 6 »i, 6 * 2 ) cos (2 (wo - w) T 12 ) + C ( 01 , 6 > 2 ) sin (2 (wo - w) T 12 )] 

dE dD, , , 

-I- / — / —- (same integrand as above ...) 

E Jdi(v) ^ 


with a{E),w,S and C as in 
delimiter than <j{E)\ 


Eqs. ( |7^ l above. 


We are now in the presence of an even narrower spectral 


X{E) = 


nh {w"^ -\- (wo ~ ‘j-')^) 


(81) 


is a (normalized) very sharply-peaked function of if = Hloq since w rin ^ ui. We can then take the 
limit piEq —>■ 00 that seems to be appropriate for the glycerol experiments and replace cr(if) by a constant, 
a{E) —>• {2hnii)~^, to get: 


AIEA {ip) 


^ d 2TT‘^fl^P* 

-Aq - XATS -- 

£{)Sr -T 0 

^<=2 dE dD , 

-^tanh 




E 


'Dr, 


E 


E^ 


2kp,T E^ - D'^ 


^nx (E) S (01,02) cos (2 (wq - w) T 12 ) 


IE, 


°° dE 
~E 


+ -TT 


rD2{D dD, . , , ^ 

/ (same integrand as above...) 

' Di(ip) 


D 


(82) 


Eor fixed and large T 12 the above is essentially the energy convolution of a function very much like the DOS 
gATs{E) times some slowly-varying functions of E and the sharply-peaked spectral delimiter x (E). The result 
of this convolution is depicted in Fig. [^and very much depends on whether Dmin < hu; or Dmin > Eo. In 
the first case, the evolution of the gATs{E) with ip (or B) gives rise to a sizable magnetic effect on the echo 
amplitude (or lEA) even in the strict-resonant limit, whilst in the second case, Dmin > the magnetic 
effect is drastically reduced and completely vanishes in the strict resonant limit, x (E) -A S{E — tiui). 

This is the explanation for the so-called isotope effect in amorphous glycerol, which actually is a mere mass- 
substitution effect. Indeed, we claim that for natural glycerol the situation Dmin > Eo applies and that for 
the deuterated samples, as for the multi-silicates, Dmin < Ejj applies instead. Then, D^in is roughly given 
by (Section 3): 


D 



moreover : Dq 


e-K'^^MVb 

V M 


(83) 


where k is the bonding constant that depends only on the chemistry and not on the isotope’s mass, M. 
Thus we immediately see that Dmm (deuterium) = Dmin (hydrogen)/-\/2 and this is sufficient to make the 
transition from case b) to case a) in Fig. 19 when a substantial mass change of the tunneling particles through 
isotopic substitution is made. The case of partial isotopic substitution requires separate considerations. A 
further consequence of mass substitution, as seen in Eq. (83), is that a larger mass M will make - the 
chemistry being unchanged - the parameters D^rnin, 

The 


much smaller and thus will give rise to a much 
slower variation of the echo’s amplitude with B, as is indeed verified in the experiments (Fig. 


14). 


highly non-uniform shape of the DOS for the ATS model thus qualitatively (and for the multi-silicates also 
quantitatively) explains all of the experimental findings. 
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Figure 19: (colour online) Convolution of the DOS gATs{E, <p) and the spectral delimiter x(^) in the integral 
giving the lEA. a) If Dmin < ^ there is a significant effect on the echo amplitude as ip increases; b) In the 
opposite case Dmin > ^ the effect is much reduced and vanishes in the strict resonant case. This is our 
qualitative explanation for the “isotope effect” [KH] . 


The last topic we discuss is our explanation for the dramatic effect for the echo’s lEA dependence on the 
waiting time ri 2 near B = 0 [59]. To explain this we take the B —>■ 0 limit of Eq. (821, when only the second 
term contributes and we get, evaluating the D-integral exactly and lumping all slowly-varying functions into 
an overall constant A*: 


/.(X 

limAIEA{(p) ^-A* / 

Jd„ 


dE 

liT 


tanh 


\2kBTj 


In 


^Omax 

^ Domin 




(84) 


xx{E) cos (2 (wo - uj) T 12 ) 


Inspection of this last expression shows, that the ATS contribution decreases quadratically as the AB phase 
(or the magnetic field) increases; hence, the disappearance of the oscillations in T 12 is slow (as experimentally 
reported jSlj) for very weak but increasing B. To understand the oscillations themselves (and their absence 
for the silicates) we set i? = 0 and redefining the overall constant we arrive at (since is the shortest decay 
time involved): 

pOO ^ — 1 

\imAIEA{0) ^-B* -^ - . 7 —r -^cos(2 (wq - w) T 12 ) (85) 

JD^i„ ^ -t (wo - w) J 

{B* being another overall constant). The remaining integral can be rewritten as an integral in the interval 
[ 0 ,cx)], which yields an exponentially decaying contribution in T 12 , plus the integral in the energy interval 
[0, hu! — Drain]- The latter is responsible for the oscillations in T 12 since: 


dx 


cosx 


x'^ -\- 2a 


Im 


e~ 6 

-Ei(—a-|-i 6 )-Ei (a- 

2 o 2a 


lb) 


( 86 ) 


where Ei(z) denotes the exponential-integral function. Rearranging Eq. (85), also adding the 2LS standard 
exponentially decaying contribution, we come to the following functional form for the total lEA of the sample: 


2 ti 2 'V / \ ^^12 

lim/E;A(0) « C 2 Lse~EAA j^CATS^^e~^) 

T 12 

f 6^ 

-I- ATATsIm I Ei (—a -I- ib) — —Ei (a -I- ib) 
2a 2a 


where C' 2 ls, Cats and Kats are appropriate constants and where: 


a — 2ri2r^^3j 


b = 2 ti2 w - 


Dr, 


(87) 


( 88 ) 


33 



























Figure 20: (colour online) The lEA from Eq. (87) plotted as a function of T 12 for parameters as in the text, 
a) The simplified theory almost reproduces (black curve) the experimental data (diamonds) for glycerol-da 
at i? = 0, nominal frequency CLi/27r=0.887 GHz and T =13.5 mK [3] (blue curve: - 2LS contribution, green 
and red curves: - ATS contributions, black curve: - sum of 2LS and ATS contributions); b) Disappearance 
of the lAE oscillations when T 0 ( 3 ) decreases or Aw increases (blue curve). From [S3]. 


We plot expression (87) in Fig. 20 as a function of ri 2 , assuming some reasonable value of r 0 ( 3 )=l.O /xs for 
the ATS spectral-diffusion time, T0(2)=1O-O for the 2LS one, and of Aw = w — HmWi = 0.6 MHz for the 
frequency offset. As is shown in Fig. 20’a) (black curve), this (over)simplified theoretical treatment almost 


reproduces the data for glycerol-d 3 at i? = 0. Fig. 20 fo) also shows (blue curve) that reducing to T 0 ( 3 )=O .3 
fjLS causes the oscillations to disappear; alternatively, this can be achieved by increasing Aw. 

This treatment also explains the findings in Fig. |18| Leaving the sample in the air causes a change in Dmin 
simply because more RERs or micro-crystals are nucleated (hence a larger Aw) by the absorbed water vapour, 
but also because of mass substitution through dissociation. The subtle role of the true value of the resonating 
frequency w (the sample sits in a reentrant resonating cavity which is loaded with a reacting dielectric and 
thus develops its own resonance frequency) is as yet not completely understood. 

Clearly, if Aw < 0 (case of glycerol-do) there can only be exponential decay (from the 2LS contribution, 
mostly), but if Aw > 0 is also too large the oscillations die out. We believe this is the case for the multi¬ 
silicates. 

Qualitatively at least, our ATS model explains all experimental observations so far. In terms of the parameters 
used for the best fits, one cannot fail to notice (see Table that for the echo experiments the values of the 
cutoffs for the combination of parameters are about one order of magnitude lower than for the other 

experiments, Cp and e, carried out inside higher temperature ranges, in fact. This can be explained by a 
mechanism whereby the number N of atomic tunneling systems within each ATS (hence within the interstices 
between the RERs or micro-crystals) picks up a temperature dependence N{T). 


7 The Magnetic Field Dependent SQUID-Magnetization 

7.1 Theory 

Having obtained qualitatively (and also quantitatively, barring the various approximations applied) good 
fits to the Cp, e and Aecho experimental data with the ATS contributions added to the 2LS’s (and when 
appropriate to the paramagnetic impurities’), we address the question of why the values of DqjS'a are so 
high. The first conclusion is that the ATSs now appear to carry considerably high magnetic moments hats 
per se. Estimating from the definition (T = 0) hats = ~ where E = ^JD'^ + is the ATS 

lower energy gap, we get for not too small fields B {hats vanishes linearly with B when B ^ 0, but saturates 
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at high enough B): 


f^j,TS-^SDo = ^(\^\SDo). 
$n (Po ^ e ) 


(89) 


Thus the very same combination \SDq of parameters appears, whilst (j)Q = hje is the electronic magnetic 
flux quantum. Using the values extracted from the Cp best fit (e.g. Table we deduce from Eq. (89) 
that (for Duran) fj,ATS ranges from about 3.8/iB to 27.1/iB. This fact alone indicates that a large group of 
correlated charged atomic particles is involved in each single ATS and that an important ATS contribution 
to the sample’s magnetization is to be expected (Fe^’*' and Fe^’*' have magnetic moment /ij = 2-\/6/rB and 
-v/Sh/Ts, respectively). 

The magnetization M of a sample containing dilute paramagnetic impurities as well as dilute magnetic-field 
sensitive ATSs is, like Cp, also given by the sum of two different contributions: 

a. Langevin’s well-known paramagnetic impurities’ contribution (Fe^“'" and Fe^“'", with nj concentration 
of one species having spin J), given by the standard expression [71] 


Mj = njg^iBJBj{z), 

where the Brillouin function Bj is defined by: 


z = 


ggBBJ \ 

kBT J 


r. ^ ^ 2J-bl , 
Bj[z) = coth 


(2J-bl) 
2J ' 


- 2J 


2J 


(90) 


(91) 


and its low-field susceptibility is the known Curie law: 


^ ^ njg gB^jJ + ^) 

B ~ ikBT 


(92) 


b. the ATS tunneling currents’ contribution, given by the following novel expression as the sum of contri¬ 
butions from ATSs of lowest gap E : 


Mats = P*nATSj^ 


dE tanh 
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dE tanh 


2kBT 
\2kBT 
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- 

TYitn 

E^ - 
E^ - 


(93) 


and which can be also re-expressed (like in the case of Cats) using y = 2 k bT following form: 


Mats = 27r P*nATskBT— 

rS 


poo 


+ 


dy tanh yin 
dy tanh y In 






(94) 


with, as before: 


Ed = and E ^2 = y/D, 

Eel,2 ^ _ _ D 


2 

min 


Dl 


' LD^ 

Omaxr 


Xcia — 


2 k bT 


^min — 


2 feBT ’ 


etc.; 


We present this expression here for the first time, also motivated by the fact that we expect a contribution 
to the measured magnetization M from the ATSs that is comparable to, or even greater than, Langevin’s 
paramagnetism of the diluted Fe impurities. The above expression follows from a straightforward application 
of standard quantum statistical mechanics, with 




nATs{— 


dB-sLS X 


uats being the ATSs’ concentration (a parameter always lumped together with P*) and with H^ls given 
by Eq. (§. The angular brackets (• • •) denote quantum, statistical and disorder averaging. 


35 



















BAS glass (B = 3 T) 



Figure 21: (colour online) The best fit of the magnetization data 
the Fe^+ and Fe^’*' impurities) and Eq. (94) (for the ATSs). 


for the BAS glass, using Eq. (90) (for 


The above formula for Mats is in fact correct for weak magnetic fields. For higher fields the correction 
implied by Eq. (20) has to be implemented and it corresponds to an improved analytic expression for the 
lowest ATS gap: 




(95) 


which needs to be used for high fields [S3]. In practice, as seen in Eq. (20), this corresponds - where 
appropriate - to the replacement of B with B^l — A(i 3 /iJ*) 2 _ ^ full derivation and the study of the B- 
and T-dependence of Mats will be presented elsewhere. Here, it suffices to say that temperature-dependent 
fitting parameters are necessary to achieve good fits for M(T) in such wide temperature range (4< T <300 
K). The underlying physical reason is that as temperature drops the RERs of the cellular structure of the 
glass (Section 2) begin to fuse together and adsorb atoms/ions from the interstices between them, growing 
at their expense. As a result the number N{T) of atomic tunnelers in each interstitial ATS decreases with 
decreeing temperature: Tg > T ^ 0 and a reasonable temperature-dependence is of the Arrhenius type: 


N{T) = No exp - 


Eq \ 

kBT] 


(96) 


where Eq is a suitable activation energy. Then, advocating coherent tunneling of the group of N{T) atomic 
tunnelers making up each interstitial ATS [15] we must use the temperature dependent parameters: 




= D. 


Eomin 

e 


Eomax 

e 


(0) 

min exp 

Eomin ~ Szs 
e 

Eomax ~ 
e 


^4 

keT 

(0) 

exp 


l(o) 


exp 


Eq \ 
ksTo J 
3Eo 
ksT 
3Eo 


+ 


3Eo 

ksEo 

3Eo 


ksT ksTo 


(97) 


and a similar one for B*, containing two parameters proportional to N{T) in Eq. (96). Here, Tq is a 


temperature corresponding to the parameters’ combinations marked with a (O)-superscript and the values we 
extract from the best fits below correspond to parameters at that temperature (typically, Tq is the average 
temperature of our Cp fits in Section 4). 


7.2 Comparison with available data 


The magnetization data [15| were best-fitted with Eq. (90) (for the Fe^'*' and Fe^"*" contributions) as well as 


with Eq. (94) (for the ATSs’), using the parameters from the Cp-fits as input. The best fit for the BAS glass 


is reported in Fig. [21] and the extracted parameters in Table dW 


36 



























Parameter 

BAS glass 

[g "1 
nFe^+ [g"^] 
P*nATS [g“^] 
Pmin [K] 
Dormn\l\S [KA^] 
Doma.\l\S [KA 2 ] 

B* [T] 

Ao [K] 

vert.offset [Am^g“^] 

1.08x10^^ 

5.01x101® 

5.74x101® 

8 . 01 x 10-2 

1.31x10® 

2.44x10® 

1.02 

0.42 

-1.04x10"® 


Table 10: Extracted parameters (from the magnetization data of [35]) for the concentration of ATS and 
Fe-impurities of the BAS glass. The vertical offset represents the residual Larmor diamagnetic contribution. 


Duran (S = 5 T) 



Figure 22: (colour online) The best fit of the magnetization data |48| of Duran, using Eq. (901 (for the Fe^’*' 
and Fe^''" impurities) and Eq. (94) (for the ATSs). 


Parameter 

Duran 

«Fe2+ [g-1] 

’^Fe3+ [g-1] 

P*nATS [g“^] 
Pmin [K] 
Dormn\^\S [KA^] 
DomacclilS [KA 2 ] 

B* [T] 

Eo [K] 

vert.offset [Am^g-i] 

3.07x1011' 
2.13x1011' 
8 . 68 x 101 ® 
5.35x10-2 
2.00 xlO® 
2.81x10® 
1.63 

0.34 

-1.97x10"® 


Table 11: Extracted parameters (from the magnetization data of |35|) for the concentration of ATS and Fe 
impurities of Duran. 
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BAS glass 

Heat Capacity fit 

Tl 

1.06x103^ g-3 = 14.23 ppm 

Tl Pe^+ 

5.00x103® g-3 = 6.69 ppm 

P*nATS 

5.19x103® g-3 

Magnetization fit 

Tl ^e2+ 

1.08x103'^ g-3 = 14.38 ppm 

Tl 

5.01x103® g-3 = 6.70 ppm 

P*nATS 

5.74x103® g-3 


Table 12: Comparison between the concentrations extracted from the two different best-fitted experimental 
data sets for the BAS glass. 


The best fit for Duran is reported in Fig. [^and the extracted parameters in Table El We recall how we 
convert the Fe-concentrations thus obtained to atomic ppm concentrations (ppma), using the formula [50] : 


nj = nj 


Na 




(98) 


where nj = Nj/M is the mass density and hj = Nj/Nat the atomic concentration of the spin-J Fe species 
(to be multiplid by 100 to get the ppm) and where: Nj is the number of Fe-ions in the sample with spin 
J the molar fraction of the i-th species, Ai the molar mass of the i-th species, Nat the total number of 
atoms in the sample and Na Avogadro’s number (6.022x10^^ mol“^) for the Fe^"*" (J=2) and Fe^+ (J=5/2) 
impurities. 


The nominal concentration of Fe^+ for the BAS glass is (using Eq. (98)): 


'^Fe%+ 


= 102 ppm 

10"® • 102 • 6.022 X 10^3 mor^ 


'^Fe%+ ~ 


80.530^ 


= 7.63 • 10^ 


(99) 


which is inadequate (as in Duran’s case below) to explain the behaviour of the heat capacity as a function of 
B, here presented (Fig. and as a function of T (studied in [T]). Table 12 summarizes the concentrations 
found from our best fits of heat capacity and magnetization data, for the BAS glass. The parameters found 
in [1] were P*nATS=6.39xlO^® g”^ and hi?e=20.44 ppm where this latter was for the Fe^"*- concentration 
only. The present study conhrms that most of the Fe-impurities in these two glasses are of the Fe^"*- type [T]. 


The nominal concentration of Fe3+ for Duran is (using Eq. (98)): 


’’^Feitn 


= 126 ppm 

_ 10"® • 126 • 6.022 X 10^3 mol”^ 
" 61 - 873 ^ 


= 1.23- lO^V^ 


( 100 ) 


which again is inadequate to explain the behaviour of the heat capacity as a function of B (see Fig. 23). 
Table summarizes the concentrations found from our best fits of heat capacity and magnetization data, 
for Duran. The parameters found in [T] were P*nATS=6.92x 10^® g“3 and fii?e=47.62 ppm where this latter 
was the Fe^+ concentration only. 

Fig. [^presents the behaviour of the different contributions to the heat capacity as a function of B for Duran; 
Cpa^+ and Cpf,i+ are given by Eq. (16), respectively with the Fe^+ and Fe3+ parameters, Cparam is the sum 
of these latter two contributions, Cats is given by Eq. (18) and the green line represents the result of the 


best fit. The dashed line corresponds to the Cp{B) one would get from the nominal concentration np^ of 126 
ppm [48] as extracted from the SQUID magnetization measurements htted with the Langevin contribution 


only (no ATS contribution). Likewise Fig. 24 presents the behaviour of the different contributions (Eq.(90) 


and Eq.( 94)) to the magnetization as a function of B, also for Duran. It can be seen that the ATS contribution 
is in both cases dominant, also (in the case of the magnetization) at the higher temperatures. We finally 
present our preliminary study of the SQUID magnetization data (also available from [?5]) for the borosilicate 
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Duran 

Heat Capacity ht 

Tl + 

3.21x10^' g-^ = 33.01 ppm 

Tl Pe3+ 

2.11x101^ g-i = 21.63 ppm 

P*nATS 

8.88x10^® g-i 

Magnetization fit 

Tl Fe'^T- 

3.07x10^' g-^ = 31.58 ppm 

Tl 

2.13x101^ g-i = 21.86 ppm 

P*nATS 

8.68x10^® g-i 


Table 13: Comparison between the concentrations extracted from the two different best-fitted experimental 
data sets for Duran. 


Duran (T = 1.23 K) 



Figure 23: (colour online) The curves represent the different terms that contribute to the heat capacity of 
Duran in our best fit of the data from |48| . The dashed curve is for Langevin’s contribution only, but with 
the nominal concentration of hp’e3+=126 ppm (no ATS). 


Duran (S = 5 T) 



Figure 24: (colour online) The curves represent the different contributions to the magnetization of Duran in 
our best fit of the data from [?5]. 
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Parameter 

BK7 

nFe^+ [g"^] 
nFe^+ [g"^] 
P*nATS [g"^] 
Pmin [K] 
Dormn\l\S [KA^] 
Domax\l\S [KA2] 

B* [T] 

Eo [K] 

vert.offset [Am^g“^] 

6.69x101^ = 0,71 ppm 

3.43x1015 = 0.36 ppm 
1.40x1016 
5.99x10-2 
8.87x10^ 
1.20x105 

1.30 

0.48 

-1.08x10-5 


Table 14: Extracted parameters (from the magnetization data of [48]) for the concentration of ATSs and Fe 
impurities of the BK7 63.530 g mol~^ |1H|)- The vertical offset represents the residual Larmor 

diamagnetic contribution. 


glass BK7, for which however no substantial magnetic effect in the heat capacity Cp has been reported [48j . 


This glass has a nominal Fe-impurity concentration of npf,3+=6 ppm |7||351|S2|, yet our best fit in Fig. 25 


with both Langevin (Eq. (j^) and ATS (Eq. (94|) contributions produces the concentrations and parameters 
given in Table 14 The best fit was carried out with knowledge of ATS parameters from our own theory m 
for the magnetic effect in the polarization-echo experiments at mK temperatures [52] . We conclude that our 
main contention is once more confirmed, in that the concentration of Fe in BK7 we extract in this way is only 
about 1.1 ppm and the bulk of the SQUID magnetization is due to the ATSs. Table [Td] reports our very first 
estimate of uatsP* for BK7. Assuming P* to be of order 1 and about the same for all glasses, we conclude 
that the concentration uats of the ATSs nesting in the RERs is very similar for all of the multi-silicate 
glasses by us studied for their remarkable magnetic effects. From the present SQUID-magnetization best 
fits we have obtained 5.74x10^® g“^ (BAS glass), 8.68x10^® g“^ (Duran) and 1.40x10^® g“^ (BK7). The 
almost negligible magnetic effect in Cp for BK7 is due, in our approach, to the low values of the cutoffs 
Domin and D^rnax for this system (these parameters appearing in the prefactor and in the integrals’ bounds 
determining the ATS contribution to Cp [H]). We find values of the combinations of parameters Derain]-S 
and DQrnax\\S again as high as those extracted from Cp, e and A^cho experimental data. We believe [H] this 
is to be ascribed to the coherent tunneling of an average number N(T) of atomic tunnelers trapped - and 
strongly interacting - in each interstitial ATSs (Fig. [^. 


BK7 glass {B = 3 T) 



Figure 25: (colour online) The best fit of the magnetization data 
and Fe^+ impurities) and Eq. (94) (for the ATSs). Data from [iS] 


for BK7, using Eq. (90) (for the Fe^+ 
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8 Nature of the Tunneling Systems and the Potential Energy Land¬ 
scape 


So far we have described how the anomalous behaviour of glasses in the low temperature regime is related to 
localized hoppings between adjacent local configurations, through qm tunneling. The STM postulates that 
in the configuration space each TS can be represented by a double well potential (DWP) and makes certain 
assumptions regarding the distribution of its asymmetry and barrier height. Although from a qualitative 
point of view this picture is certainly reasonable and appealing, in practice one has to see whether in a real 
glass-forming system such DWPs with the right characteristics are indeed present. Through experiments, 
it is very difficult to obtain microscopic information about the properties of the TSs since they do not give 
direct access to the relevant information (distribution of barrier heights, asymmetry, nature of the motion of 
the participating particles, and so on (but see the attempts in [64ll65] A. 

Thus at present the only promising approach to make progress in this field are atomistic computer simulations. 
In the past such simulations have already helped to clarify many properies of the potential energy landscape 
(PEL) of glass-forming systems and allowed to understand their relaxation dynamics [751 [721 [TU [75] . (We 
recall that the PEL is the surface of the potential energy as a function of the 3A"at coordinates of the system, 
where Nat is the number of particles. Note that this PEL does not depend on temperature, but the regions 
that the system explores with a significant probablity does.) One can therefore hope that such simulations 
will also permit to gain insight into the properties of glasses at intermediate, low, and very low temperatures. 
The question that has to be addressed in such simulations is whether or not one does indeed find in the 
PEL local DWPs that have the right features. Note that one expects that the number of local minima is 
exponentially large in Nat, be. it is a huge number and therefore it will be completely impossible to find 
all of the DWPs. Eurthermore one faces the conceptual problem of how this purely classical concept of the 
PEL can be used to describe the tunneling processes, i.e. a quantum phenomenon. The simplest approach to 
address these issues is the introduction of a reaction coordinate that parametrizes the hopping form one local 
minimum to a neighboring one and from this to calculate the effective barrier height and thus the tunneling 
parameters. Subsequently this information can be fed into theoretical approaches like those discussed in the 
previous Sections to calculate quantities like the heat capacity Cp and so on. 

Eirst of all we have to discuss the order of magnitude of the main parameters that characterize the TSs, i.e. 
(for the 2LSs) the tunneling matrix element Aq, the energy of the single well liO, and the barrier height Vq. 

The tunneling matrix element Aq is given, for example, by Eq. 


I, i.e. Ao = 


e From 


the literature one can infer the typical range of values for the fir^ two quantities for the most studied glasses 
(the multi-silicates): 


• 10 ® K < Aq < 10 ^ K for the tunneling parameter |55] . 

• 10“® K < hn < 10 K for the well energy [251 [S3] 


(it is thought that HQ should be the largest energy in the problem, hence of the order of the gap energy’s 
upper cutoff Emax)- This limits have been verified also in the context of the present ETM theory’s various 
applications and allow us to deduce the range of variation of the barrier height Vq. Figure 26 shows some 
iso-lines of Aq as a function of the barrier height Vq and the well energy hrt. The range mentioned above 
for Aq corresponds to the region delimited by the two red bold lines. This plot shows that the typical 
barrier height range goes up to a maximum of 20 K or so and therefore this value represents a discriminant 
threshold for the relevant parameters of the TSs (2LSs, specifically). In the past most simulation studies of 
the properties of the PEL have focused on the so-called “inherent structures” (ISs), i.e. the local minima 
of the PEL [75] [^ [TU [7S]. Following the pioneering work of Goldstein who was the first to point out 
that the corrugated nature of the PEL is important to understand the slow relaxation dynamics of glass¬ 
forming liquids m, researchers started to study the properties of this landscape for simple glass-forming 
systems [75] [77] [7H1 [75]. Weber and Stillinger [77] were the first to analyze the PEL by studying how 
the system evolves in this complex landscape and explores the various minima. For this they periodically 
quenched the configurations of molecular dynamics (MD) simulations to zero temperature and identified 
nearby minima as well as the connecting transition states, i.e. the saddle points (SPs). Subsequently Heuer 
and coworkers carried out detailed studies in an attempt to determine the properties of the DWPs that are 
relevant for the low temperature properties of glasses [50] [SI] ■ This was done by searching for the DWPs in 
the PEL in a wide range of the energy asymmetry A, distance between the wells d, and barrier height Vq. 
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Figure 26: Iso-lines of the tunneling parameter Ag as a function of the barrier height Vq and the single-well 
energy hU. The physically relevant range for Ag is delimited by the two bold red lines. 


After parametrizing these quantities with a polynomial function and subsequently applying an interpolation 
procedure, they inferred the properties of TSs from the collected sets of DWPs. However, their conclusions 
regarding the nature of the DWPs in terms of configurations and parameters (asymmetry,...) depend strongly 
on the manner they have scanned the PEL and as a consequence one does not have, in fact, at present a 
real good understanding on the nature of the DWPs, or more generally, of the possible atomic configurations 
that are mutually connected via a tunneling process. 

As discussed above, the barrier height of the TSs that are relevant at low T is rather small. However, most 
previous studies have found that the overwhelming majority of the barriers between two adjacent ISs is not 
that small, i.e. it is a sizable fraction of the interaction energy between the particles, thus on the order of 
100 K [78]. Hence this indicates that in fact the TSs are not related to a process in which one particle changes 
its nearest neighborhood and thus moves a good fraction of the typical nearest neighbor distance, i.e. 0.3-2 
A. So, what other options remain for excitations that give rise to DWPs that have a small barrier height 
and thus a reasonable tunnelling parameter? One possibility is that these excitation occur in the interior of 
the basin of attraction of an IS. In other words, it could be that the basin of attraction is not just a simple 
slightly deformed bowl, but is instead corrugated and characterized internally by small channels and valleys. 
These kind of local deformations in the structure of the PEL can give rise to effective DWPs and therefore 
we believe that both the TSs and the ATSs are associated not with jumps of the system from one IS to a 
neighboring one, but rather to the motion inside the basin of attraction of a given IS. The TSs could be 
therefore associated with the motion from one of these valleys to a nearby one and the tunneling can occur 
among different channels. 

To make this idea more specific we show in Fig. what is intended to be an oversimplified representation of 
the real PEL of a glass, that we computed in order to clarify the concepts; as can be noted, the immediate 
surrounding of the IS is still harmonic, but as soon as one moves away (always remaining inside the same 
basin of attraction) different valleys (ending roughly at the N, E, S and W cardinal points in the drawing) 
develop (as indicated by solid lines), separated by regions in which the potential is higher and where the 
tunneling processes could occur (a tunneling path is marked by a bold green dashed line). We have two ISs 
(ISi and IS 2 ) that are connected via a reaction path (thin dashed line) that goes through a SP. Note that 
this type of transition usually involves only a small number of degrees of freedom, i.e. on the order of 4-10 
particles, since the excitation is localized in the three dimensional space. Thus the vast majority of the other 
particles change their position only marginally in this transition event, since they basically react only in an 
elastic manner, i.e. without changing their neighborhood. As mentioned above, the typical energy of the 
SPs is too large to permit a tunneling process with a reasonable probability and thus such processes might 
be irrelevant for the properties of glasses at low temperatures. We therefore conclude that the TSs should 


42 



3 


800 

600 

400 

200 

0 

-200 


-400 


-2.5 -2 -1,5 -1 -0.5 0 

Figure 27: (colour online) Simplified view of the glass potential energy landscape (PEL) of a glass-forming 
system. We show two local minima, i.e. inherent structures (ISs), that are connected by a first-order saddle 
point (SP). Also shown are some of the valleys (bold black solid lines) that emanate from the IS and often 
lead to a SP. The bold green dotted lines represent possible transition pathways from one valley to another. 
Tunneling takes place across the “reliefs” of the PEL. Note that the pathis of these valleys do not necessarily 
lead to a saddle. 

not be thought as a restricted subset of ISs forming DWPs, but rather as small irregular valleys that exist 
inside a single IS’s basin and characterize the very bottom of the landscape. As is exemplified in Fig. the 
2LSs are tunneling processes of the N-W dashed bold green-path type, whilst the much shallower N-E and 
S-E tunneling bold green-paths are probably good examples for the more elusive ATSs. We note, however, 
that to each IS are connected on the order of Nat valleys (bold solid lines), many of which lead to one of 
the 0{Nat) SPs. These valleys are local minima within the basin of attraction of a given IS. Since these 
valleys start all at the IS, their energy difference is small, as long as they are close to the IS. Furthermore 
one can expect that also the barrier that the system has to cross when moving from one valley to another 
is not very large. (Such processes are indicated by the bold green dotted lines with arrows.) Thus such an 
inter-valley transition can approximately be described by a DWP in which the two wells have almost the 
same energy and the barrier height is not very large, i.e. the DWPs do indeed have the features that we need 
in order to describe the TSs. Thus at this stage we can speculate that the TSs are related not to the IS-IS 
transitions, how is usually done, but to valley-valley transitions inside the basin of attraction of a single IS. 
For the moment we can of course only say that the valley-valley transitions will give rise to DWPs that have 
the right semi-quantitative features. It will thus be important to check by means of numerical simulations 
to what extent these DWPs can indeed be used to describe the low-temperature anomalies of glasses on a 
quantitative basis. In addition it will be important to understand how the specific macroscopic properties 
of the glass former (composition, local structure etc.) influences the features of the PEL within the basin 
of attraction of an IS, and hence to establich a link between the macroscopic properties of the glass and its 
behaviour at low temperatures. 

From the numerical-analysis point of view, in order to understand the low temperature anomalies found 
experimentally, there is a need to perform local and systematic analysis of the bottom of the IS basin. Since 
it has been found [78| that, around the glass transition temperature, the system moves quite close to the IS 
and the majority of the saddles found is of order 1, the new approach would be the exploration of the IS’s 
surroundings starting from configurations obtained with standard MD simulations at the lowest temperature 
that can be equilibrated. The mere location of the SPs is not enough to characterize the PEL, in fact, as we 
can imagine, a given IS is associated to more than one nearby saddle and furthermore there is no guarantee 
that following the direction of the softest eigenvalue necessarily leads to the saddle, as depicted in Figure [27| 
the branches that develop rather tend to vanish in the complexity of the PEL. 
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We therefore intend to analyze the properties of simple glass-forming systems in order to systematically locate 
all of the TSs and verify whether they can represent a probe of the very local structure of the PEL and, 
therefore in this sense, contribute also to the properties of the glass transition. Such numerical simulations 
should not be too difficult and are currently on the way. 

9 Conclusions and Outlook 

In conclusion, the ATS tunneling model for the magnetic effects in multi-component glasses (the multi¬ 
silicates BAS (or AlBaSiO), Duran and BK7) and contaminated mono-component glass glycerol has been 
fully justified in terms of a (not entirely) new vision for the intermediate-range structure of real glasses. In 
this scheme the particles are organized in regions of enhanced regularity (RERs) and more mobile charged 
particles trapped in the interstices between the (on average) close-packed RERs. These are coherent (owing 
to proximity and strong Coulomb forces) atomic tunnelers that can be modeled in terms of a single quasi¬ 
particles with highly renormalized parameters. This model explains a large number of experimental data and 
facts with remarkable consistency also in terms of cross-checks like the determination of the concentration of 
trace paramagnetic impurities. 

The fact that pure amorphous silica (Spectrosil-I, for example |S2]) shows no detectable magnetic effects is a 
consequence of the extremely small size of the RERs for a-Si02, deprived of almost any nucleation centres for 
both micro-crystal and polycluster (RER) formation. These RERs will therefore trap a very small number 
N of dangling-bond ionic particles, or none at all given the covalent nature of the Si-0 bonds. Hence, no 
magnetic effects are observable in the purest, single-component glasses. 

The present theory shows that the magnetic (and compositional) effects are a mere manifestation of the 
inhomogeneous, cellular-like intermediate atomic structure of real glasses. A cellular-type structure that has 
been advocated by scientists, especially (but not only) in the ex-USSR, now for almost a century. It is not 
impossible that with this vision the TSs could become in the near future the right probes with which to 
study the structure of glasses in the laboratory. 

We have also commented on the possible atomistic computational checks for this theory and it is hoped that 
in the near future the MD simulations will provide strong evidence for the nature, here proposed, of the 
tunneling states. 
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